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ABSTRACT 

Very high-quality spectra of 24 metal-poor halo dwarfs and subgiants have been acquired with 
ESO's VLT/UVES for the purpose of determining Li isotopic abundances. The derived ID, non- 
LTE ''Li abundances from the Lii 670.8 nm line reveal a pronounced dependence on metallicity but 
with negligible scatter around this trend. Very good agreement is found between the abundances 
from the Lii 670.8 nm line and the Lii 610.4 nm line. The estimated primordial ^Li abundance is 
^Li/H = 1.1 — 1.5 • 10~^°, which is a factor of three to four lower than predicted from standard Big 
Bang nucleosynthesis with the baryon density inferred from the cosmic microwave background. Inter- 
estingly, ^Li is detected in nine of our 24 stars at the > 2(t significance level. Our observations suggest 
the existence of a ^Li plateau at the level of logeeLi « 0.8; however, taking into account predictions for 
^Li destruction during the pre-main sequence evolution tilts the plateau such that the ^Li abundances 
apparently increase with metallicity. Our most noteworthy result is the detection of ^Li in the very 
metal-poor star LP 815-43. Such a high ^Li abundance during these early Galactic epochs is very 
difficult to achieve by Galactic cosmic ray spallation and a-fusion reactions. It is concluded that 
both Li isotopes have a pre-Galactic origin. Possible ^Li production channels include proto-galactic 
shocks and late-decaying or annihilating supersymmetric particles during the era of Big Bang nucle- 
osynthesis. The presence of ^Li limits the possible degree of stellar ^Li depletion and thus sharpens 
the discrepancy with standard Big Bang nucleosynthesis. 

Subject headings: Stars: abundances - Stars: atmospheres - Stars: Population II - Galaxy: evolution 
- Cosmology: early Universe 



1. INTRODUCTION 

Lithium's two stable isotopes ~ ^Li and ^Li - con- 
tinue to pose intriguing questions for astrophysicists con- 
cerned with understanding the origins of this light ele- 
ment. These questions concern the abundance of lithium 
in cool interstellar gas and every type of star in which 
lines of neutral lithium are either detected or are poten- 
tially detectable. In this paper, we address by new ob- 
servations questions about the origin of the lithium seen 
in the atmospheres of metal-poor stars. In particular, we 
provide an adequate set of determinations of the isotopic 
abundance ratio ^Li/^Li with which to constrain some 
theoretical ideas on the origins of the rare ^Li isotope 
and of the majority ^Li isotope for which the Big Bang 
is widely identified as the initial and major origin. 

An observational link of ^Li with the Big Bang was 

'BASED ON OBSERVATIONS COLLECTED AT THE EURO- 
PEAN SOUTHERN OBSERVATORY, PARANAL, CHILE (OB- 
SERVING PROGRAMS 65.L-0131, 68.D-0091 AND 273.D-5043) 



promoted first by Spite & Spite (1982) who showed 
that the lithium abundance in the warmest metal-poor 
dwarfs was independent of metallicity for [Fe/H] < —1.5. 
Their derived abundance of lithium in the dwarfs was 
given as logeLi= 2.05 ± 0.15 on the usual scale where 
log Eh = 12.00. The constant lithium abundance defining 
what is commonly called 'the Spite plateau' suggested 
that this may be the lithium abundance in pre-Galactic 
gas provided by the Big Bang. In the succeeding years, 
many observers have reinvestigated the lithium abun- 
dance of the Spite plateau but the lithium (^Li) abun- 
dance has strayed little from the Spites' original esti- 
mate. Here, we also present a new determination of 
the plateau's lithium abundance as extrapolated to zero 
metallicity that remains faithful to the 1982 estimate. 

Early observations showed that the lithium on the 
plateau was primarily ^Li, as expected from a stan- 
dard Big Bang: Maurice et al. (1984) set an upper 
limit ^Li/^Li< 0.1 for a couple of stars. The first 
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probable detection of ^Li in a very metal-poor star was 
reported by Smith et al. (1993) for HD 84937 with 
^Li/^Li= 0.06 ± 0.03. Confirmation was provided by 
Hobbs & Thorburn (1994, 1997), Smith et al. (1998) 
and subsequently by Cayrel et al. (1999) who ob- 
tained ^Li/^Li = 0.052 ± 0.019. Smith et al. (1998) 
observed seven single stars and reported ^Li in one ad- 
ditional very metal-poor star (BD -t-26°3578, also known 
as HD 338529). Cayrel et al. provided a possible detec- 
tion of ^Li in BD +42°2667 as did Nissen et al. (2000) for 
G 271-162. Non-detections gave upper limits variously in 
the range ^Li/'^Li < 0.03 - 0.10. 

Two questions now face students of the abundances 
of ^Li and ^Li in very metal-poor stars. First, predic- 
tions of the nucleosynthesis by the Big Bang are tightly 
constrained because analysis of the anisotropies of the 
cosmic microwave background determine the one free 
parameter in a standard cosmological models that sets 
the abundances (relative to ^H) of the products ^H, 
^He, "^He, and ^Li. The free parameter is the ratio 
of baryons to photons rj = n\j/n^ or equivalently the 
fraction of the critical density 57;, provided by baryons: 
J7b = 3.652 X l{)^ri/h? where h is the Hubble constant in 
units of lOOkms^i Mpc"^ The recent WMAP-based 
analysis of the cosmic microwave background implies 
rJb/i^ = 0.0224 ± 0.0009 (Spergel et al. 2003), which 
in standard Big Bang nucleosynthesis corresponds to an 
abundance of logern^ 2.65 ± 0.10 (Coc et al. 2004; 
Cuoco et al. 2004; Cyburt 2004). This predicted abun- 
dance of ^Li is about 0.5 dex or a factor of three larger 
than the measured abundance of lithium on the Spite 
plateau. Thus, one key question has become - How does 
one bridge the 0.5 dex gap between observation and pre- 
diction? 

The second question arises from the detections of ^Li. 
Our results suggest that there may be a ®Li plateau par- 
allel to the Spite plateau for ^Li with the implication that 
the major fraction of the ^Li (and presumably some of 
the ^Li) may have been synthesised prior to the onset of 
star formation in the Galaxy. The ^Li abundance of the 
plateau exceeds by a large factor the ^Li expected from 
a standard Big Bang. Thus, our second key question is - 
How does one account for the high abundance of ^Li in 
some metal-poor stars? 

Our principal concern in undertaking the observations 
reported here was to extend high-precision measurements 
of ^Li to more metal-poor stars and to lower metallici- 
ties. The observations and the analysis techniques are de- 
scribed below in Sections 13 and 13 respectively. Section^ 
deals with the determination of stellar parameters, abun- 
dances, and macroturbulent velocities. The derivation of 
the Li isotopic abundances from the observed line shapes 
of the Li I 670.8 nm resonance line is outlined in Section 
while the analysis of the subordinate Lii 610.4 nm line 
is presented in Sectional In Sectional we take up the 
two questions just outlined: How can one reconcile the 
predicted ^Li abundance using the WMAP-based Vl\^h^ 
with the measurement of the abundance for the Spite 
plateau? How does one account for those ^Li abundances 
that greatly exceed the Big Bang prediction? 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Selection of Stars 



Programme stars were se lected from the uvby-j3 cata- 
logue of Schuster & Nissen l)1988D and from the Li abun- 
dance survey of very metal-poor stars by Ryan et al. 
fl99!l). In order to be able to study the possible evolu- 
tion of the Li isotopes in the Galactic halo, we aimed at 
a group of stars covering the metallicity range —3.0 < 
[Fe/H] < —1.0 and situated in the turnoff region of the 
Hertzsprung-Russell-diagram, where the upper convec- 
tion zone of the stars is relatively thin so that proton 
destruction of Li is minimal. Furthermore, the stars had 
to be brighter than ~ 11 in order to obtain very high 
resolution spectra with sufficiently high S/N . Finally, 
known double-lined spectroscopic binaries were avoided. 

2.2. Spectroscopic Observations 

Spectra of 23 stars were obtained with UVES (Dekker 
et al. 2000) on the ESO VLT 8.2m Kueyen telescope 
during three observing runs: July 23-25, 2000; February 
5-7, 2002, and August 28-31, 2004 (5.5 hours in service 
mode). About 1/3 of the time during the first two pe- 
riods was lost due to high humidity. The seeing varied 
from about 0.5" most of the time in July 2000 to oc- 
casionally more than 2" in February 2002. One of the 
nights in July 2000 was interrupted by "Target of Op- 
portunity" observations - later compensated by service 
observations on August 4, 2000. 

In order to obtain a spectral resolution of about R ~ 
120 000, we used image slicer #3, which for an entrance 
aperture of 1.5" x 2" produces five slices each 0.3" wide. 
Thanks to the image slicer, the efficiency of UVES was 
quite good even during periods with bad seeing. In av- 
erage conditions a S/N of about 500 per spectral bin 
(0.00268 nm at the Lii 670.8 nm line) could be obtained 
for a,V ~ 11.0 star in 2 hours. The total integration time 
was always split into at least three separate exposures so 
that cosmic rays hits on the CCD could be removed by 
comparison of the different spectra. 

The July 2000 and February 2002 spectra were taken 
in the red arm of UVES with cross disperser #3 and a 
central wavelength of 705 nm. With this non-standard 
setting, the EEV CCD covers the spectral region 600 - 
700 nm in 15 echelle orders, while the MIT CCD covers 
the region 710 - 820 nm in 11 orders. We note that the Li i 
670.8 nm line is well centered on the EEV chip and that 
the orders are separated by more than 140 pixels so that 
despite of the large width of the orders (about 70 pixels) 
due to the image slicer, there is ample interorder space 
for measuring the CCD background including scattered 
light. Furthermore, the sky spectrum could be monitored 
through an extension of the entrance slit beyond the po- 
sition of the image slicer. For further details we refer to 
Dekker et al. (g002). 

The service observations in August 2004 were carried 
out with a standard UVES setting of the red arm (cross 
disperser #:3, central wavelength 600 nm). In this set- 
ting, the EEV chip covers the 500 - 595 nm spectral re- 
gion, and the MIT from 605 to 700 nm. The primary 
purpose with this alternative setting was to check an in- 
teresting Li isotopic abundance ratio obtained for the 
very metal-poor star LP 815-43 based on the July 2000 
observations. 

Further details about the spectroscopic observations 
are given in Table ^ In addition to the 23 stars listed 
here, G 271-162 is also included in the present paper. 
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Its spectrum was obtained during the commissioning of 
UVES in October 1999 with the same setting of UVES 
as in July 2000 and February 2002 but without the use 
of the image sheer, which was not available at the time. 
Instead the star was continuously moved along the slit 
to simulate the effect of the image sheer. This s tar has 
already been analyzed by Nissen et al. l|2000() but is 
re-analyzed here for comparison purposes. 

2.3. Data Reduction 

The spectra have been reduced using standard tasks 
in the IRAF echelle package with the following steps in- 
cluded: i) bias subtraction, ii) fitting a 2-dimensional 
polynomial to the flux in the interorder regions and 
subtraction of scattered light, iii) pixel- by-pixel division 
with a normalized flat field exposed without the image 
sheer, iv) tracing and extraction of the echelle orders, v) 
corresponding extraction of a thorium calibration spec- 
trum, vi) wavelength calibration, vii) rectification of the 
echelle orders by fitting cubic spline functions with a 
wavelength scale of about 0.5 nm to the continuum, viii) 
correction for the Doppler shift of the spectral lines. 

For some stars this procedure was carried out for the 
individual sheer spectra, which were then combined to 
one spectrum by applying weights according to the ex- 
posure level of the slices. The resulting spectra were com- 
pared to spectra obtained by extraction with a long slit 
covering all five slices. The difference between the two 
sets of spectra is marginal in terms of S/N and resolu- 
tion. Hence, the simpler long slit extraction was applied 
for all stars as our final reduction procedure. Further- 
more, we note that the echelle orders are not merged, 
which means that lines near the edges of the orders, e.g. 
the oxygen triplet at 777.4 nm, can be measured in two 
consecutive orders. 

The sky background spectrum turned out to be negli- 
gible for all stars observed on July 23-25, 2000, August 4, 
2000, and February 5-6, 2002. On these dates the Moon 
was close to first or last quarter. The observations of 
LP 815-43 on August 30 - 31, 2004, however, took place 
during full Moon conditions, and the sky spectrum (cor- 
responding to the extent of the five slices) turned out to 
be about 1% of the signal in the spectrum of LP 815-43. 
Hence, in this case the sky spectrum was subtracted. 

The wavelength calibration is based on 20-40 thorium 
lines per echelle order. The rms of a third order polyno- 
mial transformation between pixel and wavelength is of 
the order of 0.12pm (1.2mA). In addition, the Th lines 
can be used to check the resolution of UVES; the intrin- 
sic width of a typical Th line is negligible compared to 
the measured width (Hobbs et al. 1999). The profiles 
of the Th lines are well approximated by Gaussians, but 
the resolving power as defined from the FWHM of the 
Th fines {R = A/AAfwhm) varies along a given spectral 
order from about 100000 to 130 000 on the EEV chip. 
As can be seen in Fig. ^ we achieved slightly higher 
resolving power during our first observing run in July 
2000 than in February 2002. The lower resolution for 
the August 2004 spectra is caused by charge diffusion in 
the MIT chip, which, as mentioned above, was used in 
that period. In any case, the large number of Th lines 
available, enable us to map the variation of the resolving 
power along an echelle order. Hence, the actual instru- 
mental broadening for a given stellar line is known and 
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Fig. 1. — Measured resolving power of UVES based on the 
FWHM of thorium hnes in the echelle order containing the Li I 
670.8 nm line. Filled circles: July 2000 (EEV CCD); Open cir- 
cles: February 2002 (EEV CCD); open squares: August 2004 (MIT 
CCD). The curves are 3rd order fits to each set of data. 



can be applied in connection with the spectrum synthe- 
sis. 

In addition to increasing the efficiency, the image sheer 
in UVES also serves to broaden the spectrum over a 
larger area of the CCD. This helps to minimize prob- 
lems in the flat-fielding, i.e., to reduce the amplitude of 
any residual fringing existing after the flat-fielding. From 
exposures of bright B-type stars during the July 2000 ob- 
servations we estimate that the residual fringing is less 
than about 0.1%. For some unknown reason, the resid- 
ual fringing turned out to be larger for the February 2002 
spectra, reaching an amplitude of about 0.3 % on a wave- 
length scale that could affect the profiles of the spectral 
lines. Hence, we decided to determine a "residual fringe 
function" from the continua of several very metal-poor 
stars observed in February 2002 taking advantage of the 
fact that these stars have only a few spectral lines, which 
are shifted to different positions in the spectra due to 
large differences in radial velocities. This residual fringe 
function is well defined and was applied to all spectra 
observed in February 2002. A similar procedure was ap- 
pfied to the spectra of LP815-43 obtained on August 30- 
31, 2004; in this case the spectrum of HD 140283 from 
August 28, 2004 was used to define the residual fringe 
function. 

As a test of the quality of the spectra and of pos- 
sible problems with the correction for residual fringing 
in particular, we compare in Fig. [21 the Li i spectra of 
HD 140283 ([Fe/H] = -2.4) obtained in July 2000 and 
Feb. 2002. As seen, there is no significant trend in the 
residuals, and the rms of the difference is 0.0016 corre- 
sponding to a S/N of about 600. This is close to what 
one would expect for the S/N values of each spectrum 
given in Table ^ A similar comparison of the spectra of 
HD 140283 from July 2000 and August 2004 does reveal 
a significant trend in the residuals, but this is mainly due 
to the lower resolution of the August 2004 spectrum. 

The S/N of the spectra in the region around the Lii 
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Fig. 2. — Comparison of the Li I region in the July 2000 spectrum 
of HD 140283 with the February 2002 spectrum (shifted by 0.05). 
The lower panel shows the difference between the two spectra. 
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Fig. 3. — Sample spectra around the Lil 670.8 nm line with the 
spectra of CD-48°2445 and CD-33°3337 shifted 0.10 and 0.20, 
respectively. Note, that the region between the Li I and the Ca I 
lines in the spectrum of CD— 33° 3337 is affected by several faint 
lines, which are not identified on the figure. 



670.8 nm line ranges from about 400 for some of the very 
metal-poor stars to more than 1000 for the relatively 
bright star HD 140283. All but two of the stars have 
S/N > 500. Representative spectra are shown in Fig. 
for a range of metallicities. 

2.4. Radial Velocities and Binary Stars 

Our conclusions regarding the lithium isotope ratio 
could be compromised if our targets are binaries. We 
have searched available databases for evidence suggest- 
ing that our program stars are not single. The extensive 
radial velocity survey of Carney et al. (1994) contains 
11 of our 23 stars. One of them, G 126-062, is reported 
as a single-lined binary; the orbital solution by Latham 
et al. (2002) gives a period of 219 days and a radial ve- 
locity amplitude oi K — 4.8 km s^^. In addition, Carney 
et al. list BD -f 03°0740 as a suspected binary. Another 
accurate survey, that of Nordstrom et al. (2004) based 
on radial velocities measured with the Coravel spectrom- 
eter, contains eight of our stars but except for G 126-062 
none of them are noted as binaries. Finally, Ryan et al. 
(1999) have nine of the most metal-poor stars in common 
with us, but again there is no evidence of radial velocity 
variations. 

To supplement these studies we have determined radial 
velocities for all our stars. The values are based on the 
measured Doppler shifts of four of the lines that are used 
in Sect. 14.51 to determine the macroturbulent velocity 
broadening. The following laboratory wavelengths from 
the NIST database are adopted: Cai 612.2217nm, Cai 
616.2173nm, Fei 623.0726 nm and Cai 643.9075 nm. In 
determining the radial velocity one should ideally correct 
for gravitional redshift and convective blueshift, but as 
the net effect is expected to be below 0.5kms~^ we have 



not made such a correction. 

The measured radial velocities (corrected to heliocen- 
tric values using the IRAF task rvcorrect) are given in 
Table The accuracy of the velocities depends on the 
stability of UVES because the thorium comparison spec- 
trum was obtained only a few times per night. Compar- 
ing the radial velocities for stars observed on different 
nights we estimate that the precision is ±0.3kms~^, i.e. 
similar or slightly better than the precision obtained in 
the radial velocity surveys mentioned above. A compari- 
son to the data in these three sources and those in the Li 
isotope study of Smith et al. (1998) gives the following 
mean differences (us - them) and rms dispersions of the 
differences (excluding G 126-062): 

Carney et aL (1994): Aw^ = +0.6 ± 0.4kms"i (10), 
Nordstrom et aL (2004): Av^ = -0.4 ± 0.6kms-i (7), 
Ryan et al. (1999): Av^ = -hO.5 ± 0.3kms-i (9), 
Smith et al. (1998): Av, = -0.4 ± 0.4kms-i (7), 

where the number within parenthesis denotes the number 
of stars in each comparison. Evidently, there are small 
systematic differences in the zero point of the radial ve- 
locities, but the small rms dispersions are as expected; 
hence, no additional evidence for binarity is revealed. In 
particular, we note that BD -|-03°0740 that was suspected 
to be a binary by Carney et al. is included in three of 
the four comparisons and shows no sign of variability. 

Out of the 23 stars in Table ^ there are only three 
(HD 213657, CD -33°3337 and CD-48°2445) for which 
our radial velocities could not be compared to other ac- 
curate sources. The spectra of these stars show no evi- 
dence of lines from a companion; the same is true for the 
other stars including G 126-062. Hence, we feel confident 
that all stellar spectra can be used for determining the 
Li isotopic ratio. 
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3. ANALYSIS: METHODS, MODELS AND INPUT DATA 
3.1. Stellar Model Atmospheres 

In the present study, we have made use of both ID LTE 
hydrostatic and time-dependent 3D hydrodynamical stel- 
lar model atmospheres. One-dimensional model atmo- 
spheres with convection treated by the mixing length 
theory for the particular stellar parameters of our sam- 
ple have been generated with the MARCS code (Gustafs- 
son et al. 1975; Asplund et al. 1997), which includes 
realistic continuous and line opacities for a wide range 
of atomic and molecular species. The radiative transfer 
in the construction of the model atmospheres is treated 
by opacity sampling. As all our stars are metal-poor, 
the chemical composition of the models have an over- 
abundance of the a-elements of [a/Fe] = -1-0.5 and a he- 
lium abundance of logenc = 10.93 (Asplund et al. 2005). 
This He abundance corresponds to a He mass fraction 
of F = 0.253, which is close to the expected value from 
standard Big Bang nucleosynthesis (e.g. Steigman 2005). 
For test purposes we have also employed theoretical ID 
Kurucz model atmospheres without convective overshoot 
as computed by Castelh et al. (1997). Those models are 
based on very much the same approximations and as- 
sumptions as our standard MARCS models and have been 
interpolated to the specific stellar parameters of our sam- 
ple from the available grid of models. 

A novel feature of our analysis compared with the pi- 
oneering works on Li isotopic abundances in metal-poor 
stars of Andersen et al. (1984), Smith et al. (1993, 1998), 
Hobbs & Thorburn (1994, 1997), Nissen et al. (1999), 
Hobbs et al. (1999) and Cayrel et al. (1999) \ is our use 
also of time-dependent 3D hydrodynamical stellar model 
atmospheres for the stars for which such 3D models with 
appropriate stellar parameters exist. The 3D model at- 
mospheres of metal-poor stars are taken from Asplund 
et al. (1999) and Asplund & Garcia Perez (2001). In 
short, the models have been constructed by solving the 
standard hydrodynamical equations of conservation of 
mass, momentum and energy coupled with a simulta- 
neous solution of the 3D radiative transfer in a repre- 
sentative volume of the stellar atmosphere and near sub- 
surface layers (e.g. Stein & Nordlund 1998; Asplund et 
al. 2000). The 3D models are based on realistic contin- 
uous and line opacities and equation-of-state computed 
for the appropriate [Fe/H] under the assumption of LTE 
and opacity binning (Nordlund 1982). At any given time 
the simulation box is sufficiently large to cover at least 
ten upflowing regions (granules). The time sequence for 
each simulation employed is at least one convective turn- 
over time divided into thirty or more individual snap- 
shots to give statistically significant results for the spa- 
tially and temporally averaged line profiles. In contrast 
to the ID model atmospheres, no mixing length theory 
enters for the calculation of the convective energy flux. 
Instead, the time-dependent hydrodynamical simulation 
automatically and self-consistently predicts the convec- 
tion efficiency as a function of depth for the given stellar 
parameters. 

^ We recognize even earlier observational studies by, for example, 
Herbig (1964), Feast (1966) and Cohen (1972) who attempted to 
derive Li isotopic abundances using both the center-of-gravity and 
profile fitting methods. These works were restricted to Galactic 
disk stars and yielded no convincing ^Li detections. 



3.2. Spectral Line Formation 

The spectral line formation is treated under the as- 
sumption of (strong) LTE making the source function 
equal to the Planck function: {t) = B^[T{t)], i.e. 
continuum scattering is treated as an absorption. The 
lines used in our study are located at wavelengths at 
which Rayleigh and electron scattering make a negligi- 
ble contribution to the source function; we have ensured 
that identical results are obtained with the code used 
herein and the code bsyn from the Uppsala spectrum 
synthesis package that treats continuum scattering prop- 
erly within the LTE framework. The radiative transfer 
is solved with a Fortran implementation of the Feautrier 
method to which is coupled an IDL routine for the 
analysis to find the optimum solution in terms of element 
abundance, wavelength shift, continuum normalization 
and extra broadening beyond the instrumental broad- 
ening (ID: macroturbulence and rotation, 3D: rotation 
only). The rotational broadening is performed through 
a disk-integration of the angle-dependent intensity pro- 
files. In the ID case, microturbulence is prescribed and 
not part of the ^^-optimization while the parameter does 
not enter the 3D calculations. The same code is able 
to handle both ID and 3D geometries with otherwise 
identical calculations, which is a great advantage when 
investigating the possible systematic errors when rely- 
ing on ID plane-parallel model atmospheres. The back- 
ground continuous opacities are taken from the MARCS 
project (Gustafsson et al. 1975; Asplund et al. 1997; 
Gustafsson et al., in preparation). We have verified that 
identical ^Li/^Li ratios and absolute Li abundances to 
within 0.01 dex would result with the independent spec- 
trum synthesis code bsyn which comes as part of the 
Uppsala package. 

The resulting velocity field from the predicted con- 
vection and waves in the 3D models introduce Doppler 
shifts, which translate to shifts and asymmetries of spec- 
tral lines. These will depend on the line formation depths 
and the atmospheric conditions (in particular tempera- 
ture, pressure and velocity) for the individual lines, mak- 
ing the detailed line shapes sensitive probes of the stellar 
atmospheres. As will become clear in Sect. 15.71 this re- 
sults in a significant improvement in the ability to repro- 
duce the detailed shapes of spectral lines (e.g. Dravins & 
Nordlund 1990; Asplund et al. 2000; Nissen et al. 2000; 
AUende Prieto et al. 2002). Since the presence of ^Li in 
the stellar atmosphere manifests itself as an asymmetry 
in the Lii 670.8 nm line in excess of the unresolved fine 
structure of the ^Li line, it is clearly important to have an 
understanding of other possible mechanisms which can 
introduce line asymmetries, such as convection. It should 
be emphasized that no microturbulence or macroturbu- 
lence is necessary for the calculation of the line formation 
in 3D. 

Our final ^Li/^Li ratios will be taken from the ID anal- 
ysis rather than from the 3D calculations. Partly, this is 
due to the simple fact that 3D models do not yet exist 
for all combinations of stellar parameters of our sam- 
ple. Another reason is that the very steep photospheric 
temperature gradients in 3D simulations of metal-poor 
stars are likely conducive to significant departures from 
LTE, at least for neutral minority species (Asplund et 
al. 1999). This has already been confirmed by detailed 
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multi-level 3D non-LTE calculations for Li (Asplund et 
al. 2003; Barklem et al. 2003) which showed significant 
over-ionization of neutral lithium but the lithium abun- 
dance from 3D non-LTE calculations was not greatly dif- 
ferent from ID non-LTE calculations. Likewise, one may 
expect similarly large 3D non-LTE effects for many of 
our calibration linos of for example Fc I and Ca i used to 
determine the intrinsic atmospheric line broadening (As- 
plund et al. 1999; Asplund 2005). Perhaps, more impor- 
tant will be the line asymmetries inherent in predictions 
for 3D model atmospheres and which will depend some- 
what on the thermodynamic conditions (LTE or non- 
LTE). Under these circumstances, the derived ^Li/^Li ra- 
tios from a 3D LTE analysis may be prone to systematic 
errors. While the above-mentioned 3D non-LTE studies 
based on only a couple of snapshots from a simulation 
arc expected to provide a good estimate of the overall Li 
line strength and thus the Li abundance, predicting the 
detailed line shape requires taking the temporal average 
of many more snapshots to obtain statistically significant 
profiles. Unfortunately, performing 3D non-LTE calcula- 
tions for long time-sequences for a stellar sample our size 
even for a relatively simple atom like Li is beyond reach 
of current computers. Therefore, the here presented 3D 
LTE calculations will mainly serve as a test of the ID 
results for the Li isotopic ratios to assess possible sys- 
tematic errors. 

3.3. Atomic Data 

The basic data for the Lii 670.8 nm resonance line 
is fortunately very well-known. The wavelengths for 
the individual hyperfine and isotope components are 
taken from the highly precise frequency measurements 
reported by Sansonetti et al. (1995), which are accu- 
rate to within 0.01pm (0.1mA). The 5/- values for the 
different ''Li and ^Li components are based on the cal- 
culations by Yan et al. (1998, sec also Yan & Drake 
1995) that give / = 0.7467871 ± 0.0000010 and / = 
0.7467527 ± 0.0000010 for the ^Li and ^Li doublets. Re- 
cent experimental values are very similar to the theoret- 
ical calculations. For example, the radiative lifetime of 
the 2p upper level (27.11 ± 0.06 ns) measured by Volz 
& Schmoranzer (1996) is very similar to the prediction 
by Yan et al. (27.117301 ± 0.000036ns). Our adopted 
wavelengths and gj-values arc the same as in Hobbs et 
al. (1999) and Smith et al. (1998). The transition prob- 
abilities and wavelengths for the isotope and hyperfine 
components for the subordinate Lii 610.4 nm were also 
taken from Sansonetti et al. (1995) and Yan et al. (1998). 
While not important for the Lii lines, we note for com- 
pleteness that we use the accurate pressure broadening 
data of Barklem et al. (2000b). 

The transition probabilities for the lines of other 
species employed in the present study are taken from 
a variety of sources, all of which are expected to be accu- 
rate and thus not introduce any severe systematic errors. 
The log g/-values for the Fe II lines arc from Bicmont et 
al. (1994) while those of the Fei linos are from O'Brian 
ct al. (1991). Both of those should be accurate to within 
±0.03 dex. The corresponding data for the Oi triplet at 
777 nm comes from Wiese et al. (1996): log 5/ = 0.369, 
0.223, and 0.002 The transition probabilities for the Cal 
lines are from Smith & Raggett (1981) and for the K i line 
{\oggf = -0.17) from Reader et al. (1980). Hyperfine 



and isotope splitting have been taken into account for the 
K I line although the effects turned out to be negligible. 
We emphasize that the exact choice of transition proba- 
bility for our calibration lines used in our determination 
of the Li isotopic abundances is of no consequence for the 
determination of the line broadening in the x^-analysis. 
The excitation potential of the lower level for the lines is 
taken from the NIST database while the radiative broad- 
ening comes from the VALD database. For neutral lines 
we rely on the pressure broadening from Barklem et al. 
(2000b) while for Fe 11 the classical Unsold theory with an 
enhancement factor of 2.0 has been adopted; the recent 
improvements in computing line broadening for Fe 11 lines 
by Barklem & Aspelund- Johansson (2005) only became 
available after our analysis was completed. We empha- 
size that as our Fell lines are all weak {W\ < 4pm), 
the exact treatment of the pressure broadening for these 
lines is of no importance for our purposes. 

4. STELLAR PARAMETERS 

4.1. Effective Temperature 

In a previous paper by Nissen et al. (2002) on oxygen 
abundances derived from the [O i] line in a subset of the 
same UVES spectra as employed here, the effective tem- 
perature was derived from the h — y and V — K colour 
indices calibrated on the infrared flux method (Alonso 
et al. 1996). The Ic relative error of Tdi was estimated 
to be about ±70 K with the largest contribution com- 
ing from the uncertainty in the interstellar reddening 
A{E{b — y)) ~ ±0.01. Since it is of great importance 
for our discussion of the total Li abundance in the stars 
to obtain a higher precision of Tog, we derive here the ef- 
fective temperature from the profile of the Ha line, which 
is very sensitive to temperature and practically indepen- 
dent of interstellar reddening. 

The Ha line is well centered in one of the echelle or- 
ders on the EEV CCD. The main problem in getting an 
accurate profile is to determine the echelle blaze func- 
tion of the Ha order. The method of fitting cubic spline 
functions to the continuum, which works well for narrow, 
well isolated lines like the Li i line, cannnot be applied to 
a broad line like Ha. Hence, we made a special reduc- 
tion of the Ha echelle order taking advantage of the fact 
that the blaze function varies smoothly and slowly with 
echelle order. This allows us to determine the blaze func- 
tion of the Ha order by interpolating in the blaze func- 
tions of the two adjacent orders. A similar technique has 
been successfully applied by Barklem et al. (2002) and 
Korn (2002) for various echelle spectrographs, including 
UVES. 

The Ha profile is affected by a few metallic lines and 
several telluric H2O lines, but due to the high resolution 
of the spectra these narrow absorption lines can easily 
be removed by interpolating the Ha profile over the dis- 
turbing lines. 

The observed Ha profiles have been compared to theo- 
retical profiles computed for MARCS model atmospheres 
using the Uppsala LTE spectrum synthesis code bsyn. 
Stark broadening is based on calculations of Stehle & 
Hutcheon (1999) with the model microfield method, and 
self-broadening of hydrogen from Barklem et al. (2000a) 
is included. The profiles are computed for models having 
the same gravity and metallicity as the star but with a 
range of effective temperatures embracing the value es- 
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Fig. 4.— Observed H« profile for BD +03°0740 (dots) compared 
to theoretical profiles calculated for atmospheric models with the 
gravity and metallicity of the star (log 3= 4.04, [Fe/H] = —2.65) 
and three values of the effective temperature: T^ff = 6170 K, 6270 K 
(best fit), and 6370 K. The spectral bands used in the definition of 
the Ha index is marked. 
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timated from the colour of the star. An improved T^g is 
then determined from x^-fitting to the observed profile 
(see Fig. 0J). The region within ±0.6 nm from the center 
is not included in the fit because the core of Ha is af- 
fected by deviations from LTE (Przybilla & Butler 2004; 
Barklem, in preparation). Furthermore, we note that 
the continuum setting is considered as a free parameter 
in the fit, because the true continuum is ill-defined due 
to the limited width of the Ha echelle order. As an al- 
ternative to the x^-fitting we have also determined Toff 
from an Ha index defined as the ratio between the flux 
in the two regions marked by 'C in Fig. ^to the flux in 
the regions marked by 'L'. Comparison with model at- 
mosphere calculations of the index then yields Tcs- As 
with the x^-fitting, the derived Toff does not depend on 
the continuum normalization. 

The two methods to determine Tcff give consistent re- 
sults within ±20 K. We also note that the Toff derived for 
a given star from spectra observed on different nights is 
extremely stable with variations being less than ±15 K. 
The calculated Ha profile is relatively insensitive to er- 
rors in the gravity and the metallicity of the models. Al- 
together, we estimate that relative temperatures for our 
sample of metal-poor turnoff stars are determined with 
a Icr precision of about ±30 K. As discussed in Sect. 14.31 
the good agreement between Fe abundances derived from 
Fe I and Fe II lines, respectively, suggests that this small 
error estimate is realistic. 

The absolute temperature scale is, however, more un- 
certain than ±30 K, as it depends on the theory of hy- 
drogen line broadening among other things. Had instead 
the older Ali & Griem (1966) theory of resonance broad- 
ening been applied instead of the more complete self- 
broadening theory of Barklem et al. (2000a), the Tcfi de- 
rived would be raised by about 100 K. In addition, there 
are uncertainties in the absolute Toff values of at least 
100 K arising from the current generation of ID model 
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Fig. 5. — Difference in Tcff as estimated from Ho and from (b — y) 
and {V — K) photometry calibrated on the infrared fiux method as 
a function of T^fj {upper panel) and [Fe/H] (lower panel). 



atmospheres, which may not capture the effects of con- 
vection in a realistic manner. Detailed statistical equi- 
librium calculations for Hi based on available e~+Hi 
and Hi+Hi collision cross-sections and ID model atmo- 
spheres reveal no departures from LTE for the wings of 
Ha while the line cores are much better described (Przy- 
billa & Butler 2004; Barklem, in preparation). For com- 
pleteness, we note that the adopted He abundance has 
no effect on the Ha-based Teg: a 0.2 dex lower He abun- 
dance, as perhaps induced by elemental diffusion in the 
star, would only increase Toff by about 8K. 

The temperatures derived from Ha are given in Table 
12 When compared to the effective temperatures in Nis- 
sen et al. (2002, 2004) determined from b-y a,ndV - K 
calibrated on the infrared flux method, good agreement 
for the majority of the stars is found: the mean difference 
is AToff = — 34±95K. A few stars show large differences: 
the Ha-based Toff for G 126-062 is 170 K higher while for 
two stars (G 020-024 and BD -13°3442) the photometric 
temperatures are higher by about 200 K. The reason for 
the latter is probably an overestimate of the reddening 
for the two stars in Nissen et al. (2002, 2004). OveraU 
the agreement is quite satisfactory, and there are no ob- 
vious trends in the temperature differences with Tcff or 
[Fe/H], as seen in Fig. We are thus unable to find any 
evidence from Ha for the very hot Tcg-scale of Melendez 
& Rami'rcz (2004) and Ramirez & Melendez (2005a,b) 
based on their reapplication of the infrared flux method. 



8 



Asplund et al. 



4.2. Surface Gravity and Absolute Magnitude 

The surface gravity is determined from the fundamen- 
tal relation 

logA^log^+41og;^+ (1) 

OA{Mboi - Mhoi,Q) 

where Ai is the mass of the star and Mtoi the absolute 
bolometric magnitude. 

The absolute visual magnitude My is derived from the 
Hipparcos parallax (ESA 1997), if available, and also 
from Stromgren photometry using the method and cali- 
brations in Nissen et al. (2002, 2004). The only excep- 
tion from the procedure described by Nissen et al. is that 
we used the spectroscopic Ha index (Sect. 14.1(1 instead 
of the photometric H/3 index to derive the interstellar 
reddening excess. This provides a more precise estimate 
of E{b — y), and consequently the limit for performing a 
correction of the Stromgren photometry for interstellar 
absorption and reddening was set at E{b — y) > 0.010 
instead o{E{b-y) > 0.015 in Nissen et al. (2002, 2004). 

If the parallax is available with an accuracy Att/tt < 
0.3, the mean value of Mv,phot and Mv,par is adopted 
for My, otherwise My^phot- The bolometric correction 
is taken from Alonso et al. (1995), and the mass of a 
star is derived by comparing its position in the My- 
logTcff diagram with evolutionary tracks calculated by 
VandenBerg et al. (2000). 

As the determination of log g depends on Tcff directly 
and [Fe/H] indirectly the determination of the stellar pa- 
rameters has to be iterated. The final values are given 
in Table El together with values for [Fe/H] and [0/H] 
(Sect. 14.31 and Sect. 14. 4() and the microturbulence ^turb 
(Sect. I4.5|l . The internal error of the derived value of 
logg is dominated by the error of My, which we have 
estimated to be AMy = ±0.20. This corresponds to 
Alog5= iO.OSdex. While the VandenBerg et al. (2000) 
models do not include the effects of diffusion, this trans- 
lates to the surface gravity possibly being overestimated 
by at most 0.04 dex (e.g. Profhtt & Michaud 1991). We 
note that this effect is much too small to impact any of 
our conclusions in what follows. 

The upper panel in Fig. El shows the My- log Tcff dia- 
gram for stars with [Fe/H]< —1.7 compared to two sets 
of evolutionary tracks from VandenBerg et al. (2000) 
corresponding to [Fe/H]= —2.31 (the mean metallicity 
of the stars) and —2.01, respectively, and an a-element 
enhancement of [a/Fe] — +0.3. Except HD 19445 (a 
main sequence star) and HD 140283 (a subgiant), the 
stars are close to the turnoff. Some of the scatter in 
Tcff around the turnoff can be attributed to differences 
in [Fe/H] with the most metal-poor stars being warmest. 
The lower panel of Fig. shows the location of the eight 
stars with [Fe/H]> —1.7 in the HR-diagram. They have 
a mean metallicity of [Fe/H]= —1.33. Again all stars 
except perhaps HD 106038 are classified as turnoff stars. 
As discussed in Sect. 17.11 HD 106038 has a peculiar high 
Li abundance. 

4.3. Metallicity 

Iron abundances have been determined as described in 
detail by Nissen et al. (2002) from 13 Fell lines, six of 
which are in our spectra. The remaining seven lines lie 
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Fig. 6. — HR-diagram for stars with [Fe/H]< —1.7 (Upper panel) 
and for stars with [Fe/H]> —1.7 {Lower panel). Filled circles refer 
to stars with a > 2o- detection of ^Li, and open circles to stars 
with no clear detection of ®Li. The displayed evolutionary tracks 
are from VandenBerg et al. (2000). 



in the spectral region 490 - 550 nm and were measured 
from UVES spectra with i? ~ 55 000 and S/N ~ 300 
obtained with the primary aim of determining beryl- 
lium abundances (Primas et al., in preparation). For 
the very metal-poor stars BD-f09°2190, BD-f03°0740 
and BD -13°3442 only two of the Fell hues (519.76 and 
523.46 nm) could be measured reliably. For these stars, 
we added the somewhat stronger Fell 492.32 nm line to 
improve the precision of [Fe/H]. For the most metal-poor 
star in our sample, CD — 33°1173, we do not have a spec- 
trum in the 490 - 550 nm region. Hence, [Fe/H] could 
not be determined from Fe ii lines, but was estimated 
from Fei lines taking into account the offset of 0.08 dex 
between [Fe/H] from Fell and Fel lines discussed below. 
The Fe il lines were also measured in the solar flux spec- 
trum (Kurucz et. al. 1984), allowing differential iron 
abundances [Fe/H] to be determined without knowledge 
of the exact oscillator strengths of the lines. We note that 
for an adopted solar iron abundance of logeFc= 7.50, 
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the ^/-values derived from the solar spectrum agree well 
with those of Biemont et al. (1991) using the MARCS 
solar model atmosphere. 

The Fe abundances for the Sun and our sample of stars 
were determined from the observed Fe ii line strengths 
using the Uppsala eqwidth program together with ID 
MARCS model atmospheres having the appropriate stel- 
lar parameters. As T^s, logg and [Fe/H] are interlinked, 
the determination of [Fe/H] is an iterative procedure, 
with the final values listed in Table |21 The microturbu- 
lence ^turb used in these calculations was adjusted such 
that the derived Fe abundances show no trend with line 
strength. For [Fe/H] < —2.0, the employed Fell lines are 
weak and, hence, the derived Fe abundances are practi- 
cally independent of ^turb- For those stars, a microtur- 
bulence of 1.5 km s^^ has been adopted. 

The error in [Fe/H] induced by the uncertainties in Teg 
is small: A[Fe/H] < ±0.02 dex. Since Fell is the domi- 
nant ionization stage in the line formation region of these 
stars, the derived Fe ii based abundance is dependent on 
logg: an error in log 5 of ±0.1 translates to an error in 
[Fe/H] of ±0.035. Adding in quadrature the error indi- 
cated by the line-to-line scatter with those arising from 
the uncertainties in Tgg and log 5 implies a statistical er- 
ror of ±0.05 dex in [Fe/H] . The systematic error in [Fe/H] 
may be somewhat larger. 

In addition, we have measured Fe abundances based 
on up to 11 weak Fei lines using ^/-values from O 'Brian 
et al. (1991). As shown in Fig. [3 the Fel lines system- 
atically give [Fe/H] values slightly lower than those from 
Fell lines. The average difference, [Feii/H]- [Fei/H], is 
0.08 dex with an rms scatter of only 0.041 dex. The ex- 
pected scatter arising from the adopted errors, ATcff = 
±30K and Alogg= ±0.08dex, is 0.046 dex, i.e., not 
significantly different from the observed scatter. This 
suggests that the estimated precision of T^ff and log 5 is 
realistic. 

The difference in the derived Fe abundances from Fe i 
and Fell lines could reflect departures from LTE, most 
likely for Fei (see discussion in Asplund 2005) or er- 
roneous stellar parameters (e.g., too low by about 
100 K). The small abundance differences and the lack of 
a trend in [Fe ii/H]-[Fe i/H] with either T^s or [Fe/H] is 
a strong argument against pronounced departures from 
LTE for ID atmospheres, as claimed for example by 
Thevenin & Idiart (1999). 

Our [Fe/H] values are systematically higher by about 
0.2 dex than those of Ryan et al. (1999) for the 11 stars in 
common. Ryan et al. adopted literature values of [Fe/H] 
mainly based on Fe i lines, with a supporting role played 
by a calibration of medium resolution spectra. It is likely 
that this 0.2 dex difference is largely attributable to the 
discrepancy between Fei and Fell lines mentioned above, 
the lower Teff-scale of Ryan et al. and/or differences in 
the adopted g/-values. Given the remaining uncertain- 
ties surrounding Fe I lines in terms of possible departures 
from LTE, we believe it is wise to rely on Fell lines for 
the estimate of [Fe/H]. 

4.4. Oxygen abundance 

Oxygen may be a preferred alternative to iron in track- 
ing the evolution of the lithium isotopes in the early 
Galaxy. Oxygen is synthcsised primarily in the most 
massive stars and distributed when those stars explode 
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Fig. 7. — Differences in [Fc/H] estimated from Fell and Fei lines 
as a function of T(,ff {upper panel) and [Fe/H] {lower panel). While 
there is an offset of about 0.1 dex, there are no obvious trends with 
stellar parameters. This difference may signal differential non-LTE 
effects for Fe I relative to the Sun or slightly too low adopted Tcff 
by about 100 K. 



as Type II supernovae; ^Li may also be made in these su- 
pernovae by the z^-process (Woosley et al. 1990). Also, 
when spallation is an effective mode of lithium synthesis, 
spallation of oxygen nuclei is the dominant mechanism. 
Iron is also a product of Type II supernovae but in con- 
trast to oxygen its yield is dependent on the mass cut 
between the stellar remnant (neutron star or black hole) 
and the ejecta and iron is also a product of Type la su- 
pernovae. The latter factor is not a major concern here 
as our stars have a metallicity below that at which Type 
la supernovae are generally considered to begin their con- 
tributions. 

The O abundances for our program stars have been de- 
rived from the O i triplet at 777 nm. Unfortunately, the 
alternative abundance indicator, the [Oi] line at 630 nm, 
is undetectable for [Fe/H] < —2.5, even with the very 
high S/N of our spectra. We have, therefore, opted for 
using the same diagnostic for the entire sample. The O i 
triplet typically yields somewhat higher O abundances 
than the forbidden [Oi] line at low metallicity (see dis- 
cussion in Asplund 2005 and references therein). 

The Oi 777.1, 777.4 and 777.5nm lines are located 
towards the edges of two consecutive echelle orders from 
which the equivalent widths were measured and averaged 
(weighted by their respective S/N). The results for the 
July 2000 spectra have already been presented in Nis- 
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sen et al. (2002) to which we now add nine stars from 
February 2002. From the equivalent widths the ID LTE 
O abundances were derived using the stellar parameters 
listed in Table |2J it should be noted that the O abun- 
dances for the July 2000 targets differ slightly from those 
given in Nissen et al. (2002) due to our use of Ha-based 
Tgff rather than the photometric values employed in our 
previous study. We have also computed ID non-LTE 
abundance corrections for the O i triplet lines in an iden- 
tical fashion to the calculations presented in Nissen et al. 
(2002). These non-LTE corrections are always negative 
and typically —0.1... — 0.2 dex but reach — 0.3dex at the 
lowest [Fe/H]. Our final adopted O abundances given in 
Table 121 include these non-LTE corrections. These abun- 
dances are referenced relative to the same solar O abun- 
dance from a ID solar model atmosphere as adopted in 
Nissen et al. (2002), i.e. logeo = 8.74, which is 0.08 dex 
higher than the current best estimate based on 3D hy- 
drodynamical model atmospheres (Asplund et al. 2004, 
2005). 

Nissen et al. (2005a) have discussed in some detail 
the trend and scatter of [0/Fe] as a function of [Fe/H] 
for our sample. Here we just note that [0/Fe] increases 
from about 0.4 dex at [Fe/H] = -1.0 to 0.6 dex at [Fe/H] 
= —3.0. The scatter in [0/Fc] around a linear fit to the 
data ([0/Fe]= -0.11 [Fc/H] + 0.31) is 0.08 dex with a 
large contribution coming from a relatively low [0/Fe] 
in HD 160617, a star that is known to belong to the rare 
class of nitrogen-rich subdwarfs (Bessell & Norris 1982). 

4.5. Spectral Line Broadening 

In order to extract the lithium isotopic ratio from the 
profile of Li i 670.8 nm line, it is paramount to have an ac- 
curate understanding of the intrinsic spectral line broad- 
ening. This broadening comes in various forms all of 
which must be taken into account: thermal, pressure, 
radiative, rotational and Doppler broadening. Thermal 
broadening can be directly computed from the known 
temperature structure of the stellar atmosphere while 
the pressure and radiative broadening are, in principle 
at least, known atomic properties specific for each tran- 
sitions. 

It is well-known that the observed spectral lines are 
always broader than predicted solely from the thermal 
and atomic broadening. All stars are expected to ro- 
tate, which broadens the line profiles but does not mod- 
ify the line strength. As our observations are unable 
to disentangle readily rotational broadening from other 
atmospheric broadening agents, we adopt in all cases 
i;rotsini= 0.5kms~^ given the great age of these metal- 
poor halo stars (for comparison, the Sun has Wrot = 
1.8kms~^ at the equator). We have confirmed that our 
derived results are immune to this particular choice of 
Urotsini. In addition, atmospheric velocity fields from 
convection, wave motion, turbulence and the like intro- 
duce Doppler shifts which affect the line profiles. In anal- 
yses based on ID hydrostatic model atmospheres, one 
normally introduces two free parameters: microturbu- 
lence ^turb and macroturbulence Cmacro, acting on small 
and large spatial scales, respectively. Following standard 
practice we determine ^turb by ensuring that the derived 
abundances from Fe I and Fe ii lines are independent of 
line strength while Cmacro is estimated from the observed 
line shapes. The form of Cmacro is assumed to be Gaus- 



sian. This choice was made mainly for convenience as 
the effects of Cmacro Can then be easily combined with 
the line broadening due to the finite instrumental resolu- 
tion which is also described by a Gaussian distribution. 
Wc have performed various tests ensuring that our re- 
sults in terms of Li isotopic abundances are not affected 
by the use of a Gaussian Cmacro instead of, for example, 
a radial-tangential version (e.g. Gray 1992). 

The lines used for determining Cmacro were K i 
769. 9nm, Cai 612. 2nm, Cai 616. 2nm, Cai 643.9nm, 
Fei 623.0 nm. Fell 624.7 nm and Fell 645.6 nm. In some 
cases we also added the Fei 610.2 nm and Fei 613.6 nm 
lines when there was a lack of calibration lines of suit- 
able strength compared with the Lii line. For the Au- 
gust 2004 observations of LP 815-43 and HD 140283, ad- 
ditional lines could be used since the spectra extended 
further to the blue, enabling the Fell 501.8 nm, Fei 
523.3 nm, Fei 532.4 nm, Fei 539.7 nm, Fei 538.3 nm, Fei 
543. 4nm, Fei 542. 4nm, Mgl 552. 8nm, Fei 561. 2nm, 
Cai 558.9 nm and Fei 558.7 nm lines to be covered in 
exchange for the Ki 769.9 nm line. The Ki line was 
broader than other lines, such as Fe i 623.0 nm, of compa- 
rable strength in most stars. This may signal significant 
departures from LTE which have indeed been predicted 
for this resonance line for metal-poor stars (Ivanova & 
Shimanskii 2000; Takeda et al. 2002). Good agreement 
was found, however, when the K i line had an equivalent 
width of about 2 pm, i.e., similar to that of the Lii line 
itself. This is consistent with the finding of Smith et al. 
(2001) in the case of HD 84937. We, therefore, decided 
not to include this line in the final estimates of Cmacro and 
used it only as a consistency check on the mean value and 
scatter implied by the other lines. 

When computing the average broadening parameters 
from the calibration lines, we have considered only lines 
with strengths 0.9 < Wx < 7.5 pm. The calibration lines 
included in the computation of the mean broadening are 
of similar line strength to the Lii 670.8 nm line, which 
should minimize possible systematic differences. This 
point is important since we have noticed that the de- 
rived Cmacro is dependent on the strength of the line, as 
illustrated in Fig. (Sjfor the case of HD 106038 (a Li- rich 
star with a significantly stronger Lii 670.8 nm line than 
in the typical star of our sample: W\ = 6.32 pm rather 
than « 2 pm). In this particular instance, Cmacro was 
computed as a mean of the four lines with W\ = 6±1 pm, 
which is expected to give a good estimate of the intrinsic 
broadening for the Li i line. 

5. LI ISOTOPIC ABUNDANCES FROM THE Lii 670.8 NM 
RESONANCE LINE 

5.1. Li Isotopic Abundances 

The presence of ^Li alongwith ''Li in the stellar at- 
mosphere introduces additional width and asymmetry 
of the line profile^. As outlined in Sect. 14.51 we first 
determine Cmacro from a number of Ki, Cai, Fei and 

^ We will not consider here the center-of-gravity method which 
is based on the wavelength shift introduced by the addition of ^Li. 
The nearby Cai 671.7 nm line has served as a reference line. As 
discussed by Smith et al. (1998) and Cayrel et al. (1999), the 
presence of differential line shifts arising from convection (see, for 
example, Reddy et al. 2002), and the remaining uncertainties in 
the laboratory wavelength of the Ca I line render this method very 
prone to systematic errors. 
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Fig. 8. — Estimated macroturbulence ijmacro (assumed to be de- 
scribed by a Gaussian) from different lines of varying strengths in 
HD 106038; in addition the Unes were broadened by rotation (as- 
sumed to be t;rotsini= 0.5kms~^) and instrumental broadening 
(measured from Th lines). Note the strong correlation between 
Cmacro and equivalent width, which emphasizes the importance of 
having calibration Unes of similar strengths as the Li I 670.8 nm 
line when estimating the intrinsic line broadening. The equivalent 
width for the Li line in HD 106038 is denoted by the vertical dotted 
line. 



Fe II lines. This is achieved by a x^-analysis of the pre- 
dicted and observed line profiles allowing the element 
abundance, wavelength zero-point and continuum nor- 
malization to vary besides Cmacro- In most cases, Cmacro 

could be well determined with a standard deviation of 
ACmacro= ±0.1 — 0.2kms~^. Larger uncertainties of 
±0.3 — 0.6kms~^ resulted for four stars: G 126-062, 
CD-33°1173, HD 59392 and BD-h09°2190. The esti- 
mated broadening parameters are listed in Table|31 while 
Figs. ISh 121 show the resulting profiles for one of these cal- 
ibration lines (Cai 612.2 nm) for a few stars. 

Equipped with the estimate of the macroturbulence 
Cmacro for each star, the adopted Wrotsini (0.5kms~^) 
and the measured instrumental broadening at 670.8 nm 
from Th lines, synthetic profiles of the Lii 670.8 nm line 
for different Li isotopic ratios were generated for the pro- 
gram stars using the stellar parameters listed in Table HJ 
Calculations were done for ^Li/^Li between —0.10 and 
0.10 by step of 0.01. Negative ^Li/^Li ratios, of course, 
do not exist in reality but were included in order to ob- 
tain estimates as accurate as possible for stars with very 
small or vanishing ^Li content. These were produced by 
allowing the ^Li components of the Li I 670.8 nm line to 
be in emission rather than in absorption as is the case 
for positive ^Li/^Li ratios. At the low metallicity of the 
program stars, one does not need to worry about poten- 
tial blends from various atomic and molecular lines that 
can perturb the Li feature at solar metallicity (Nissen et 
al. 1999; Israehan et al. 2001, 2003; Reddy et al. 2002) 
and we thus computed the theoretical profile taking into 
account only the various components of the Li line. Test 
calculations revealed that this was fully justifiable for the 
target stars. 

The comparison between the theoretical and observed 
profiles was quantified using a x^-analysis similar to that 
used in most recent investigations of ^Li/^Li in late- 
type stars (e.g. Smith et al. 1998; Nissen et al. 1999, 
2000; Israelian et al. 2001, 2003; Reddy et al. 2002; 
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Fig. 9. — Upper panel: Observed Cal 612.2 nm profile (rhombs) 
in HD 140283 together with the best fit theoretical ID LTE pro- 
file (solid line). Also shown are the predictions with the macro- 
turbulence (^macro ±0.5kms-l. The Gaussian FWHM is here the 
combination of Cmacro and the instrumental broadening; the pro- 
files have also been convolved with a stellar rotation velocity of 
'UYot = 0.5kms~^. Middle panel: Computed Lil 670. 8nm profile 
for three different values for ®Li/^Li (0.00, 0.05 and 0.10) shown to- 
gether with the observed line (rhombs). In each case, the Li abun- 
dance, wavelength shift and continuum level have been allowed to 
vary. The location of the various '^Li and "Li components are indi- 
cated. Lower panel: Resulting Ax^ for the combined Gaussian of 
the mean Cmacro and instrumental broadening; also shown are the 
cases for fmacro iA(^macro. In this star, no significant detection of 
^Li can be claimed. 
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Fig. 10.— Same as Fig. |51 but for G 013-009. In this star, 
the x^-analysis indicate a > 2(t detection of ®Li at the level of 
6Li/''Li= 0.05. 
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Fig. 11.— Same as Fig. ^hut for G 020-024. In this star, the x^- 
analysis indicate a 4(7 detection of ^Li at the level of '^Li/^Li = 0.07. 



Aoki et al. 2004). The was computed according to 

= X) ^ /^^: where 0\ and Si denote the ob- 
served and synthetic flux at wavelength point i, respec- 
tively, and (7 = (S/N) ^ as estimated in three nearby 
continuum windows. Typically 26 wavelength points 



were considered across the Li line to estimate x^- For 
each ^Li/''Li, the total Li abundance, the wavelength ze- 
ropoint of the observed spectrum and the continuum level 
were allowed to vary in order to optimize the flt and thus 
minimize x^. The most probable value for ^Li/^Li was 
obtained by cubic spline interpolation of the resulting 
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Fig. 12.— Same as Fig. [glbut for LP 815-43. In this star, the y^- 
analysis indicate a 2a detection of ®Li at the level of ^Li/^Li = 0.05. 
The observed spectrum shown here is the one from August 2004. 



values, from which the minimum could be de- 

termined. The IfT, la and 3a confidence limits of the 
determinations correspond to Ax^ — l<d ~ Xmin — 
and 9, respectively. These uncertainties are determined 
by the noise of the spectra, which as mentioned above are 
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Fig. 13.— Derived ^Li/^'Li as a function of T^^ ( Upper panel) 
and [Fe/H] [Middle panel). The stars considered to have a signifi- 
cant detection (> 2a-) of ''Li are shown as solid circles while non- 
detections are plotted as open circles. Note that the non-detections 
do not cluster around ^Li/^Li= 0.00 but rather at si 0.01. This 
is also illustrated as a histogram {Lower panel). The observed dis- 
tribution is clearly not representative of a Gaussian centered at 
6Li/'^Li= 0.00. 



better than S/N > 500 for all but two stars that have 
400 and 470, respectively. The reduced x? is in all cases 
close to one with the variations between stars as expected 
from statistical fluctuations. In addition, the x^-analysis 
was repeated for each star with different adopted macro- 
turbulence parameters, Cmacro ± A^macro, where ACmacro 

was given by the standard deviation from the various Ca 
and Fe lines. Examples of the observed and predicted 
Li I line profiles and the corresponding x^ values can be 
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found in Figs. l^l2lfor a few stars. The resulting ^Li/^Li 
ratios and "^Li abundances are listed in Table |2I 

The errors in the derived ^Li/^Li ratios are quantified 
by taking the uncertainties from the finite S/N as well as 
the adopted Cmacro and stellar parameters in quadrature. 
In essentially all cases, the first source dominates the er- 
ror budget in spite of the extremely high S/N achieved in 
these observations, but for a couple of stars with poorly 
determined Cmacro such as CD— 33°1173 the second un- 
certainty drives the final error. The errors arising from 
the uncertainties in the stellar parameters are in all cases 
very minor, as discussed in detail in Sect. 15.51 We note 
that our error estimates are slightly more conservative 
than previous studies of the Li isotopic ratios in metal- 
poor stars, which, for example, often have relied on only 
one or two calibration lines to determine the line broad- 
ening. 

As seen in Table |31 nine of the 24 stars fulfill our crite- 
ria of a > 2cr result to be considered a significant ^Li de- 
tection: HD 102200, HD 106038, HD 160617, G 013-009, 
G 020-024, CD-30°18140, CD -33°3337, CD -48°2445 
and LP 815-43. This should be contrasted with the previ- 
ous single confirmed detection of ®Li in HD 84937 (Smith 
et al. 1993, 1998; Hobbs & Thorburn 1997; Cayrel et 
al. 1999). The stars with detected ^Li are all very near 
the turn-off for their respective metallicities (Fig. |SJ), in 
accordance with the suggestion of Smith et al. (1993, 
1998). In addition, we find a large number of border- 
line cases at the 1 — 2cr level that we do not consider 
significant enough to claim a positive detection. In fact, 
there is a tendency for the non-detections to cluster at 
^Li/'^Liw 0.01, as illustrated in Fig. US This could be 
taken as a warning sign that our analysis contains an 
unidentified systematic error that bias the derived iso- 
topic ratio towards too high values. We have therefore 
carried out tests in an attempt to identify such a prob- 
lem, which are described in detail in the following subsec- 
tions. In summary, no serious problem that could affect 
our analysis has been uncovered: no differences are found 
when taking into account possible ISM absorption in the 
Li line (Sect. 15. 2|) . adopting an alternative x^-analysis 
(Sect. I5.4|l . changing the stellar parameters within rea- 
sonable bounds (Sect. 15. 5|) . employing a different suite of 
ID model atmospheres (Sect. I5.6|l or using 3D hydrody- 
namical model atmospheres (Sect. I5.7|l . Furthermore, we 
independently obtain the same results as some previous 
investigations when relying on their observations (Sect. 
I5.8|l . Finally, our confidence in the derived results is 
boosted by the analysis of HD 19445, which was included 
in our program to provide a consistency check (see, also, 
Smith et al. 1993). Since this star is a relatively cool 
(Toff = 5980) metal-poor dwarf {logg= 4.42), standard 
models of stellar evolution predict that ^Li should be 
thoroughly destroyed on the star attaining the main se- 
quence (see discussion in Sect. I7.3|l . Indeed, an analysis 
identical to that for all of our other stars indicate no 
presence of ^Li: ^Li/'^Li = 0.002 ± 0.008. 

We conclude that the Li isotopic ratios presented in 
Table |31 are reliable. We claim significant ^Li detections 
only in the above-mentioned nine stars, but we can not 
rule out the possibility that the majority of turn-off halo 
stars in fact have ^Li/^Li « 0.01 (Fig. I13|l . To confirm an 
isotopic ratio as low as 0.01 for an individual star is ex- 
tremely challenging (c.f. our analysis of HD 140283 with 



S/N w 1000), but confirmation should be attemptable 
in at least a statistical sense for a sample like ours but 
with even higher S/N spectra and additional calibration 
lines to nail down the line broadening accurately. 

The Li abundances derived from the profile fitting of 
the Li I 670.8 nm line have been obtained under the as- 
sumption of LTE. We have applied the non-LTE abun- 
dance corrections computed by Carlsson et al. (1994) 
interpolated to the stellar parameters of our sample to 
the LTE Li abundances. These non-LTE calculations are 
particularly appropriate as they have been performed us- 
ing similar MARCS model atmospheres as employed here. 
The resulting non-LTE ^Li abundances are given in Table 
O The non-LTE abundance corrections are in all cases 
very small (AlogeLi< 0.02 dex) and generally negative. 
The scatter in Li abundances is not significantly modified 
by the application of these non-LTE corrections. 

5.2. Is the High ^Li/Lt in LP 815-43 Real? 

Due its low [Fe/H], LP 815-43 will feature prominently 
in the discussions in Sect. [7| It is natural to inquire 
whether its high ^Li abundance determined from our 
X^-analysis is real. We note that both the spectrum of 
LP 815-43 obtained in July 2000 and that acquired in Au- 
gust 2004 imply a detection of ^Li at the 2a level or more. 
The former indicated ^Li/'^Li = 0.078±0.033 based on an 
observed spectrum with S/N = 540 but with only three 
Cai lines to estimate the intrinsic line broadening. This 
result prompted us to request additional observing time 
for this star to improve the S/N and include more cal- 
ibration lines. The combined spectrum from two nights 
in August 2004 has S/N = 630 as well as covering as 
many as 10 calibration lines of the right strength which 
importantly beat down the uncertainty in the macrotur- 
bulence to 0.19kms~^. The resulting isotopic ratio is 
slightly less but still implies a significant detection by 
our criteria: ^Li/^Li^ 0.046 ± 0.022. The result from 
the July 2000 spectrum is fully consistent with this im- 
proved value taking the uncertainties into account. The 
individual spectra obtained on August 30 and 31, 2004, 
have S/N = 420 and 490, respectively. Both of these 
imply °Li/^Li« 0.03 — 0.04 but each allows only a Icr- 
detection due to the reduced S/N. Finally, we note that 
had we adopted the better determined macroturbulence 
value from the August 2004 spectrum for the analysis of 
the July 2000 spectrum, the derived ®Li/^Li would have 
increased to 0.089 ± 0.027, i.e. a 3cr-detection. 

As described in Sect. 12.31 there may be residual fring- 
ing in the reduced spectrum at the level of < 0.1% of 
the continuum level with a typical wavelength scale of 
0.1-0.2nm. If it is as high as 0.1% and one of those de- 
pressions falls on top of the Li line, it could mean that the 
^Li/^Li has been overestimated. An ad-hoc adjustment 
of the local continuum level by this amount would change 
^Li/^Li for LP815-43 to only a la detection and thus not 
significant by our standards (of course if instead of a de- 
pression there is a peak in the residual fringing it would 
increase ®Li/^Li by the corresponding amount). This 
illustrates clearly the challenges posed by this method 
to determine the Li isotopic ratio in metal-poor turn-off 
stars. An enlargement of the observed August 2004 spec- 
trum is shown in Fig. ^1 While there may be hints of a 
residual fringing at some wavelengths, there arc no traces 
thereof in other parts of the spectrum near the Li line. 
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Fig. 14. — Enlargement of the August 2004 spectrum of LP 815- 
43 around the Li I 670.8 nm hne. Note that the spike at 670.655 is 
a cosmic ray hit. Also shown with a dashed line is the predicted 
contribution of ^Li to the Li I 670.78 nm line for ^Li/'^Li= +0.05. 



We have tried to quantify the amount of fringing by 
computing the equivalent width directly integrated over 
a 0.08 nm window (the width of the Li line) along the 
echelle order with the Li line. The wavelength variation 
of these equivalent width measurements is slightly more 
pronounced than expected from pure random noise with 
the same S/N, suggesting the presence of fringing with 
a quasi-regular wavelength dependence. Adding a sinu- 
soidal component with an amplitude of 0.05% and wave- 
length scale of about 0.17 nm to such pure random noise 
produces a similar equivalent width distribution to that 
observed in August 2004, as shown in Fig. fringing 
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Fig. 15. — Histogram of the measured equivalent widths over a 
0.08 nm window in the continuum near the Li I 670.8 nm line for the 
August 2004 spectrum of LP 815-43 (solid line). Also shown is the 
same distribution for pure random noise with the same S/N (630) 
as the observed spectrum (dotted line). The dashed line shows the 
distribution of equivalent widths for the random noise spectrum 
with an added sinusoidal component of amplitude 0.05% (sec text 
for details). The similarity with the observed distribution suggests 
the presence of residual fringing at the 0.05% level for this star 
while fringing at the 0.1% level can be ruled out (not shown). 



as large as 0.1% can clearly be ruled out. A similar exer- 
cise for the July 2000 spectrum suggests residual fringing 
with an amplitude of 0.08%. It can also be noted that 
the expected ^Li contribution to the equivalent width 
(0.1pm for ^Li/^Li= -1-0.046) is very unlikely to be pro- 
duced purely by chance as the observed equivalent width 
distribution of the nearby continuum have a typical un- 
certainty of 0.03 pm, consistent with our claimed 2a de- 
tection for this star. We stress that the presence of fring- 
ing at the 0.05% level does not automatically imply that 
our derived ^Li/'^Li for this star is overestimated as it 
also requires that a fringing depression rather than a 
peak falls on top of the Li line, something we are not 
able to evaluate the likelihood of. The fact that inde- 
pendent observations acquired in 2000 and 2004 both 
yield a significant presence of ^Li in the star increases 
our confidence that our ^Li detection in this star is real. 
We nevertheless strongly encourage further observations 
of this star to confirm our conclusions. 

A further concern in the analysis of LP 815-43 is that 
lithium interstellar (ISM) absorption line may mimic 
the presence of ®Li considering that the Ki 769.8 nm 
is clearly distorted by a blue-shifted (— 2.07kms~^) ISM 
hne. Welty & Hobbs (2001) found in their survey of ISM 
absorption that N{Ul)/N{Ki) = 0.0054. Together with 
the difference in (^/-values of the two resonance lines, 
the ISM Li I profile should be about 100 times weaker 
than the corresponding ISM Ki hne. The ISM Ki line 
has a line depth of about 8% and an equivalent width 
of 0.97 pm. The ISM Lii line should, therefore, have 
an equivalent width of less than 0.01pm (0.1mA). It IS, 
therefore, unlikely that an ISM lithium line will have a 
significant infiuence on the derived ®Li/''Li. 

To quantify the above estimates, we first subtracted 
from the K i 769.8 nm profile the expected stellar K i line 
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as predicted using an K abundance similar to the other 
four stars with similar [Fe/H] and Cmacro and Urad from 
the Cai and Fei lines. We then fitted a Gaussian to this 
ISM Ki profile, which was subsequently converted to a 
corresponding ISM Li i profile while taking into account 
the hyperfine structure of the Li line and the difference 
in overall line strength. The difference in line width of 
the Gaussians due to different thermal broadening of Li 
and K was also included. After subtraction of this blue- 
shifted ISM profile, we again performed a x^-analysis 
of the remaining stellar Li i line. As expected from the 
above estimates, the resulting ^Li/^Li was indistinguish- 
able from the previous value. We conclude that our de- 
rived high ^Li/^Li is not an artifact of ISM Li absorption. 

5.3. Estimating Uncertainties Using Monte Carlo 
Calculations 

As an independent means of estimating the uncertain- 
ties attached to our Li isotopic abundance analysis, we 
have carried out extensive Monte Carlo simulations for 
HD 140283. A theoretical profile of the Lii 670.8 nm 
line with an adopted ^Li/''Li= 0.02 was treated as an 
observed line, to which noise corresponding to various 
S/N values (400, 500, 600, 700, 800, 900 and 1000) were 
added. All line profiles were convolved with a Gaus- 
sian of width 5kms~^ and projected stellar rotational 
velocity of 0.5kms~^ prior to introducing the effects of 
finite S/N. These artificial observed lines had the same 
wavelength sampling as our UVES spectra. We then 
performed a x^-analysis of these profiles in an identical 
fashion as for our real observations with the exception 
that the line broadening parameters were assumed to be 
known beforehand. For each S/N, 100 realisations of the 
noisy synthetic profile were generated and exposed to the 
X^-analysis. The scatter in the estimated ^Li/^Li gives 
an estimate of the accuracy with which the isotopic ratio 
can be determined for a given S/N. To these uncertain- 
ties should still be added the possible errors arising from 
uncertainties in the line broadening, continuum place- 
ment and stellar parameters. 

The estimated ^Li/^Li deviated from the input value 
with a standard deviation of 0.010, 0.007, 0.005 and 0.004 
for S/N = 400, 600, 800 and 1000, respectively. These es- 
timates are consistent with the results using the true 
observed profile for this star when ignoring the uncer- 
tainty stemming from the estimated errors in the adopted 
macroturbulence and stellar parameters. We conclude 
that the error estimates in Sect. 15. H are indeed realistic. 

5.4. Effects of Choice of a Analysis 

In their analysis of the Li isotopic ratio in HD 84937, 
Cayrel et al. (1999) adopted a slightly different technique 
than the one employed by us. Rather than determin- 
ing the intrinsic line broadening due to macroturbulence 
and rotation from lines of other species, they introduced 
this as an additional free parameter in the x^-analysis of 
the Li I 670.8 nm line. Only after the analysis did they 
confirm that the derived broadening parameter was con- 
sistent with values estimated from the Cal 616.2 and 
671. 7nm lines. The ^Li/^Li ratio Cayrel et al. obtained 
for this star (0.052 ± 0.019 at la) is in agreement with 
the value (0.06 ± 0.03 at la) of Smith et al. (1998) from 
an independent observational data set. Smith et al. used 
the same approach as us but with fewer calibration lines. 



We have investigated how our conclusions would be 
modified by employing Cayrel et al.'s approach for our 
full sample. Our results are reasonably unchanged by 
adopting this alternative approach. Only in six of our 
24 stars is the derived ^Li/^Li ratio changed by more 
than ±0.02. Of our nine detections at the > 2a level, 
seven (HD 102200, HD 106038, HD 160617, G 020-024, 
CD-30°18140, CD-33°3373, LP 815-43) are still de- 
tections under the same criterion. While G 013-009 
and CD —48° 2445 would no longer be significant de- 
tections, three other stars would in this framework be 
classified as > 2a detections: HD 213657, HD 298986 
and BD -|-03°0740, which otherwise were only la results. 
BD+03°0740 ([Fe/H]= -2.65) is noteworthy given its 
low [Fe/H]. 

We believe that Cayrel et al.'s method is inferior to 
ours as it does not make full use of the information con- 
tent in the observed spectra. Due to the limited wave- 
length coverage of their spectra, Cayrel et al. had access 
to only two calibration lines (Cai 616.2 and 671.7 nm) 
which in HD 84937 are of quite different strength com- 
pared with the Li line. Under those circumstances it may 
be justified to allow also the broadening parameters to 
vary in the x^-analysis of the Lii 670.8 nm profile when 
the S/N is very high as in their spectrum. One of the 
great advantages of the VLT/UVES is the wide wave- 
length coverage that provides several lines for estimating 
the intrinsic line broadening. We also note that the un- 
certainties in the estimated ^Li/^Li ratios are in general 
larger when Cayrel et al. 's method is used in place of our 
standard analyis presented in Sect. 15.11 in some cases sig- 
nificantly so (e.g. BD-13°3442: 6Li/7Li= 0.001 ± 0.028 
and 0.042 ±0.050, respectively), while in no case was the 
uncertainty reduced. We conclude that the method of 
Cayrel et al. is less appropriate than the one employed 
in the present study. Nevertheless, our results would 
have been quite similar had we adopted this approach. 

5.5. Effects of Choice of Stellar Parameters 

As described in Sect. 01 the estimated uncertain- 
ties in the adopted stellar parameters in absolute terms 
are ATcff^ ±100 K, Alogg= ±0.2, A[Fe/H]= ±0.2 and 
A^turb= ±0.2kms~^. We have repeated our Li iso- 
topic abundance determinations for all stars using MARCS 
model atmospheres with modified stellar parameters by 
deriving new line broadening parameters and then re- 
peating the x^-analysis of the Lii 670.8 nm profile. The 
resulting ^Li/^Li ratios when changing the stellar param- 
eters by +100 K, +0.2 dex and +0.5kms~^, respectively, 
are given m Table H Clearly, the ^Li/'^U ratios are ba- 
sically immune to the uncertainties at this level in the 
stellar parameters: the slightly different line broadening 
parameters derived from the calibration lines also apply 
to the Lii line, leaving ^Li/^Li essentially intact. Very 
rarely is A^Li/^Li as large as 0.003 with the difference 
being < 0.001 in the majority of cases. This finding is in 
line with those reported previously by Smith et al. (1993, 
1998). 

While the effects of the uncertainties in log 5, [Fe/H] 
and ^turb are negligible for the derived "^Li and ^Li abso- 
lute abundances (AlogeLi< ±0.01 dex for our estimated 
stellar parameter uncertainties), this is not true when 
changing Tcff. Our calculations reveal that logeLi in- 
creases with about 0.073 dex for ATcs= +100 K, neces- 
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sarily quite similar to the estimate by Ryan et al. (1999). 

5.6. Effects of Choice of Suite of ID Model Atmospheres 

The tendency for our derived ^Li/''Li ratios to be 
slightly but positively offset from zero, may raise concern 
that our analysis suffer from an unidentified systematic 
error. One possibility could be our use of MARCS model 
atmospheres, since Nissen et al. (1999) found a smaller 
isotopic ratio (0.033 instead of 0.041) for the disk star 
HD 68284 when switching from a MARCS to a Kurucz 
model atmosphere with identical stellar parameters. To 
investigate whether the choice of suite of ID model atmo- 
spheres has any effect on the derived Li isotopic abun- 
dances, we have repeated the analysis for all stars us- 
ing Kurucz model atmospheres without convective over- 
shoot and for the same stellar parameters, as computed 
by CastelH et al. (1997). The resuhing ^Li/'^Li ratios 
are presented in Table0]together with the corresponding 
values with MARCS model atmospheres. The changes to 
the isotopic ratios resulting from substitution of the Ku- 
rucz for the MARCS models are, as anticipated, negligible: 
A^Li/^Li< ±0.001 and in no case larger than 0.003. Fur- 
thermore, the small differences can be both positive and 
negative. 

Similarly, the effect of the Kurucz models on the ab- 
solute Li abundances is marginal"^. The average differ- 
ence compared with the MARCS based results is Alog 6^= 
0.015 ± 0.018 dex with the maximum difference being 
0.048 for HD 102200. 

For completeness, we note that the above test is based 
on the assumption that the stellar parameters of the 
model atmospheres need no adjustments when switch- 
ing to a different suite of model atmospheres. Given the 
almost identical behaviour of our Lii, Ki, Cai, Fei and 
Fell lines in the two cases, we expect insignificant dif- 
ferences for both the photometric calibrations or the Ha 
line profile fitting. In summary, our conclusions about 
the lithium isotopic ratio are immune and the lithium 
abundance insensitive to the particular brand of theoret- 
ical ID, hydrostatic, LTE model atmospheres employed 
in the analysis. 

5.7. Effects of 3D Hydrodynamical Model Atmospheres 

In addition to the tests described above, we have 
also performed a similar x^-analysis using the time- 
dependent 3D hydrodynamical model atmospheres de- 
scribed in Sect. 13. II under the assumption of LTE. Since 
such 3D models still only exist for relatively few stellar 
parameter combinations, only some of our target stars 
could be exposed to a 3D analysis. The 3D models as- 
signed to the various stars were based on similarity in 
stellar parameters; we note therefore that the here pre- 
sented 3D results are obtained with somewhat different 
stellar parameters than adopted for the ID calculations 

^ Differences in absolute Li abundances approaching O.fdex 
would have been found had we resorted to the convective over- 
shoot models of Kurucz (1993) instead of using the models with 
this option switched off (Ryan et al. 1996). These overshoot models 
are, however, significantly less successful in reproducing observed 
stellar flux distributions and their shallower temperature gradients 
find no support in our detailed 3D hydrodynamical simulations of 
stellar atmospheres with surface convection. We conclude that the 
overshoot models are less appropriate for abundance analyses of 
halo stars. 



described in Sect. 15.11 As explained in Sect. 13.21 the 
main advantages with a 3D approach is that neither 
micro- nor macroturbulence enter the analysis and that 
the intrinsic line asymmetries arising from convective 
motions follow directly from the 3D line formation calcu- 
lations. The projected stellar rotation velocity Vrotsini 
must however be determined from other lines. Thus, in 
the x^-analysis, v^ot sin i replaces Cmacro as a free variable; 
the Gaussian instrumental broadening measured for each 
line is known from Th lines. 

The results for the stars for which a 3D analysis has 
been carried out are given in Table while Fig. El 
and Fig. El show examples of the achieved profile fits 
for the Ca 1 616.21 lines for a few of the stars. It is 
noteworthy that the 3D analysis yields a substantial im- 
provement in the agreement with the observed profiles 
of Ca and Fe lines, as also previously reported by Nis- 
sen et al. (2000). In almost all cases, the resulting 
for the calibration lines are smaller than for the cor- 
responding ID cases; the only exceptions are for the 
K 1 769.89, Fell 624.76 and/or Fell 645.64 nm lines for a 
few of the stars. The derived Vrot sin i values are be- 
tween 0.1 and 3.8kms~^. The typical uncertainties in 
Wfot s'mi are 0.3 km s^^ but in some cases can be as large 
as « lkms~^. Very little is known about the observed 
rotation properties of metal-poor halo stars but due to 
their high age one would expect them to rotate slowly 
(the Sun has Vrot = 1.8kms~^). It is therefore quite 
possible that the estimated rotational velocities towards 
the higher end are overestimated. For example. Smith 
et al. (1998) suggest Vrot sini < 2 — 3kms^^ for a typ- 
ical metal-poor turn-off star. It can, however, not be 
ruled out at present that some of these stars actually 
have Drotsiniw 3kms^^. Part of the explanation for 
some of the high estimated Uiot sin i values is that the 
adopted stellar parameters for the 3D models are inap- 
propriate for some of the stars: if a too low Toff or too 
high log g are used then the convective motions and thus 
the Doppler line broadening will be too small, requiring 
large v^ot sin i to compensate when fitting the observed 
spectral lines. Furthermore, we remind the reader that 
LTE has been assumed for these 3D line formation cal- 
culations while the steep temperature gradients in low 
metallicity 3D models can be expected to result in sig- 
nificant non-LTE effects (Asplund et al. 1999; Asplund 
2005). We suspect that this may be the explanation for 
the slight correlation between v^t sini and [Fe/H] in our 
sample, although we note that iirotsini is particularly 
uncertain for the two stars with [Fe/H] « —1 (HD3567 
and G 075-031) because of the quite strong calibration 
lines employed. We therefore urge these Wrotsini values 
not to be taken too literally at this stage. In fact, it is 
encouraging that the resulting rotational velocities are 
at least reasonable given that there are no broadening 
mechanisms like micro- and macroturbulence as in ID 
that can be tuned. 

Interestingly, in general the 3D LTE calculations do 
not improve the agreement with the Li 1 670.78 nm line 
as judged by the x^-statistics. For more than half of 
the stars, the 3D analysis yield slightly higher minimum 
X^-values than in ID. The most likely explanation in our 
opinion for this contrasting behaviour compared with the 
Ca and Fe lines is that LTE is a particularly poor as- 
sumption for the Li i line in metal-poor stars; pronounced 
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Fig. 16.— Cal 616.21 nm line in HD 140283 (upper panel) and 
HD 160617 [lower panel). The solid lines correspond to the best 
fitting 3D profiles from a x^analysis while the dashed lines are for 
the ID profiles; the rhombs are the observed profiles. Also shown 
are the residuals between observed and predicted line profiles. 



over-ionization of Li i is predicted from 3D non-LTE cal- 
culations due the steep temperature gradients of the 3D 
model atmospheres (Asplund et al. 2003; Barklem et al. 
2003). 

An inspection of Table [S] reveal that the ID ^Li/^Li ra- 
tios appear to be similar or lower than the ratios derived 
from 3D LTE calculations. Only in three stars are the 3D 
results noticeably smaller than for the corresponding ID 
case. The smaller ^Li/^Li ratios in HD3567 and GG75- 
031 are likely a reflection of the relatively high [Fe/H] of 
these two stars, which renders the calibration lines used 
to determine Vrot sin i quite strong (equivalent widths up 
to 12 pm); in the case of the Sun, the impressive agree- 
ment between predicted and observed line profiles found 
for weak lines break down for such strong lines (Asplund 
et al. 2000). The difference between 3D and ID for 
LP 815-43 is of particular interest given the prominent 
status this star will have when later discussing the im- 
plications of our ^Li survey. In 3D, the derived isotopic 
ratio is ^Li/^Li= 0.036 ±0.021 to be compared with the 
ID value of ^Li/'^Li^ 0.046 ± 0.022 (both results are for 
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Fig. 17.— Same as Fig. QUbut for G 020-024 {upper panel) and 
G 075-031 {lower panel). 



the August 2004 observations). Taking these 3D LTE 
calculations at face value, LP 815-43 would no longer be 
considered a significant ^Li detection, although it would 
not be far from making the 2a criteria. 

Considering the expected severe 3D non-LTE effects on 
the Li I line in metal-poor stars that are not accounted 
for here, the ID results are to be preferred. Indeed, one 
would expect also the calibration lines used to estimate 
Wrot sin i to be vulnerable to departures from LTE, in par- 
ticular minority species like Cai and Fei (Asplund 2005). 
Carrying out detailed 3D non-LTE calculations for both 
the calibration lines and the Li i line for the necessary 
many simulation snapshots is unfortunately extremely 
computationally demanding. We will therefore have to 
postpone investigating this remaining issue to a future 
date. The preliminary 3D LTE analysis carried out here, 
however, has not uncovered any significant systematic 
errors that would suggest that our ID-based Li isotopic 
ratios are overestimated in general. 

5.8. Comparison with Previous Studies 

A number of our targets have been investigated for 
their Li isotopic abundances in previous studies. Our 
non-detection of ^Li is consistent with the upper lim- 
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Fig. 18.— Li I 670.78 nm line in G 020-024 {upper panel) 
and LP815-43 (lower panel) based on a 3D LTE analysis. The 
best fitting line profiles as judged by a x'^a^ii^'lysis for ®Li/''Li = 
+0.00, +0.05 and +0.10 are shown. 



its reported for HD 19445, HD 140283 and BD +03°0740 
by Smith et al. (1993, 1998) and Hobbs & Thor- 
burn (1994, 1997). Our derived value for G 271-162 
{^U/'^U= 0.019 ±0.012) does not constitute a detection 
of ^Li but is very close to Nissen et al.'s (2000) result of 
^Li/^Li = 0.02±0.01 based on the same VLT/UVES com- 
missioning spectra, in particular as wc have here adopted 
a slightly more conservative estimate of the uncertainties. 

More noteworthy is the result that we do not con- 
firm the ^Li detection for HD 338529 (also known as 
BD+26°3578) by Smith ct al. (1998) who found 
SLi/^Li= 0.05 ±0.03. Our value of 0.010 ±0.013 is based 
on a higher quality spectrum {S/N ~ 520 vs 420) but the 
main difference can be traced to the adopted line broad- 
ening parameters. Smith ct al. combined the effects of 
instrumental broadening and macroturbulence by con- 
volving the synthetic spectra with a Gaussian profile, 
but used only the Cai 671.78 nm line to determine the 
width of this Gaussian; however, this line is quite weak 
in HD 338529 having a line depth of only 3% at i? = 
110,000. As a result. Smith et al. could not constrain 
the macroturbulent velocity parameter for this star very 



well (this can be seen in their Figure 6 for HD 338529, 
where the a-s a function of macroturbulent velocity is 
very shallow). Our higher broadening parameter is sig- 
nificantly better determined as it is based on four Cai 
and Fei lines, all in very good agreement. We would also 
have found a significant ^Li/''Li ratio had we adopted the 
same macroturbulence as Smith ct al. (1998). In addi- 
tion, for two of our stars with significant ^Li detections, 
Smith et al. could only obtain upper limits: < 0.06 and 
< 0.05 (2(7) for HD 102200 and HD 160617, respectively. 
Again, the main difference appears to be the adopted 
macroturbulence parameter, which in our case is smaller 
than their value. Also, our spectra for these stars have 
significantly higher S/N. The differences in *^'Li/'^Li be- 
tween us and Smith et al. for HD 338529 certainly serve 
as a warning when trying to interpret the results in terms 
of Li isotopic abundances for this extremely challenging 
method. It is clearly paramount to have exceptionally 
high S/N as well as well-determined intrinsic broadening 
parameters based on a multitude of spectral lines with 
similar strength as the Li I 670.8 nm line. While we have 
tried to accomplish both of these factors and believe that 
we have set a new standard in this respect, we urge the 
reader to exercise some caution when using our quoted 
Li isotopic abundances. 

Deliyannis & Ryan (2000) reported a firm detection of 
^Li in the subgiant HD 140283: ^U/'^U = 0.040 ± 0.015. 
No details of this analysis have as yet appeared in the 
literature. We are, therefore, not able to identify reasons 
for the discrepancy with our non-detection of ^Li/ ""Li = 
0.008 ± 0.006. Our result for HD 140283 is, however, 
consistent with the upper limit of ^Li/''Li< 0.018 esti- 
mated recently by Aoki et al. (2004) based on a Sub- 
aru/HDS spectrum of the quality of our VLT/UVES 
spectrum. Wako Aoki has kindly made available their 
Subaru spectrum of HD 140283 for an independent anal- 
ysis by us. As we did not have access to their measured 
instrumental profile, we performed the analysis adopting 
a combined Gaussian for the instrumental and macro- 
turbulence broadening, a case also considered by Aoki et 
al. Our result of ^U/^U= 0.003 is in perfect agreement 
with their corresponding result of 0.002 when using their 
choice of calibration lines. We note that when we per- 
form an analysis in the same way as for our own sample 
(i.e. using the same calibration lines and taking into ac- 
count uncertainties in line broadening, S/N and stellar 
parameters), we find ^Li/^Li= 0.006 ± 0.004, which is 
in excellent agreement with the UVES result. When ex- 
panding the set of calibration lines to include in total 14 
lines between 500 and 650 nm from the Subaru spectrum, 
the 1(7 error doubles: '5Li/'^Li= 0.007 ±0.009. The exact 
results and uncertainties thus depend (slightly) on the 
particular choice of lines to estimate the line broadening, 
which must be borne in mind when comparing different 
^Li/''Li-analyses. 

6. LI ABUNDANCES FROM THE Li I 610.4 NM 
SUBORDINATE LINE 

6.1. Analysis 

The exceptionally high quality of our UVES spectra 
allows Li abundances to be derived from the subordi- 
nate Li I 610.36 nm line, which in these halo stars has 
an equivalent width of just 0.1 — 0.2 pm (1 — 2mA). We 
have performed a similar x^-analysis of a wavelength re- 
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Fig. 19. — Observed and theoretical Lil 610.4 nm profiles for 
HD 19445 (upper panel) and HD 140283 (lower panel). The four 
theoretical ID LTE profiles are computed with no Li, log el; = 
2.0, 2.2 and 2.4, respectively, with the stellar parameters listed in 
Table |^ For each star, the residuals between the synthesis and 
observations are also shown. 
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Fig. 20.— Same as for Fig. Q^but for LP 815-43 (upper panel) 
and BD+09°2190 (lower panel). No significant detection of the 
Lil 610.4 nm line is made in the latter star according to a x^-test, 
although there is a possible feature at the correct wavelength that 
would correspond to log^j^i^ 2.0. 



gion around 610.4 nm, taking into account the hyperfine 
and isotope splitting for the Li hne and neighboring hues 
of Fei, Fell and Cai. The adopted Fe and Ca abun- 
dances were those imphed by the other considered Fe 
and Ca hnes. While no attempt was undertaken to opti- 
mize the agreement between observations and predictions 
for these Fe and Ca lines by adjusting the gf-values or 
elemental abundances, we emphasize that this had no 
impact on the derived Li abundances as the neighboring 
lines are quite well described and sufficiently far away 
from the Li feature not to cause any problems. Fig. [T^ 
and Fig. I^OIcompare the synthesized profiles for different 
Li abundances with the observed spectra for a few stars. 

The derived Li abundances from the subordinate Li i 
610.4 nm line are listed in TableOl The line is unambigu- 
ously detected in 22 of the 24 program stars with the two 
exceptions being the very metal-poor stars BD -|-03°0740 
and BD +09°2190. Also for those two stars there is a tan- 
talizing suggestion of a feature at the correct wavelength, 
although the x^-analysis does not indicate a significant 
detection. As a consequence, we only give 3cr upper lim- 



its to the Li abundance in these two stars. 

The derived ID LTE Li abundances from the Lil 
610.4 nm line have been corrected for departures from 
LTE using the detailed ID statistical equilibrium calcu- 
lations of Carlsson et al. (1994). The non-LTE abun- 
dance corrections for the subordinate line are positive 
and slightly larger than for the resonance lines: typically 
AlogELi— -1-0.05 dex. The maximum non-LTE effect oc- 
curs for the metal-poor subgiant HD 140283 (+0.09 dex). 

6.2. Comparison with Results from the Li Resonance 

Line 

Figure 1211 compares the derived Li abundances from 
the weak subordinate Lil 610.4 nm line with those from 
the stronger Lil 670.8 nm resonance line. The agree- 
ment is very encouraging with a mean difference of only 
0.05 ± 0.05 dex for the 22 stars with measured Li abun- 
dances. For all practical purposes, the two lines return 
the same lithium abundance for a given star. The dif- 
ferences show no obvious trends with metallicity or ef- 
fective temperature. The upper limits for the Li abun- 
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Fig. 21. — Differences in the ID non-LTE abundances from the 
Li I 610.4 nm and 670.8 nm lines as a function of T(,ft (upper panel) 
and [Fc/H] [lower panel). 



dances from the subordinate line in BD -|-03°0740 and 
BD -|-09°2190 are very close to the abundances from the 
resonance line. It is unlikely that the explanation for 
this minor abundance difference can be found in erro- 
neous ^/-values for either of the two lines. One (remote) 
possibility is that the adopted T^s values are underesti- 
mated; a « 150 K higher Toff-scale would erase the mean 
difference. Another explanation could be inappropriate 
non-LTE abundance corrections for the ID models. Use 
of 3D models with non-LTE line formation might erase 
the differences (Asplund et al. 2003). A study aimed 
at investigating the 3D non-LTE line formation for a 
wider range of stellar parameters as well as revisiting the 
ID non-LTE abundance corrections in light of recent im- 
provements in coUisional data (Belyaev & Barklem 2003; 
Barklem et al. 2003) would be worthwhile. 

The Li i 610.4 nm line was first detected in the halo star 
HD 140283 by Bonifacio & Molaro (1998). They also 
found very good agreement in the abundances inferred 
from the two Li lines. Likewise, Ford et al. (2002) de- 
rived the same abundances from the 610.4 and 670.8 nm 
lines for this star using a very high quality spectrum ob- 
tained with Subaru/HDS (the same spectrum as used by 
Aoki et al. (2004) to derive an upper limit to ^Li/''Li 
in this star as discussed above). However, for the other 
seven stars for which Ford et al. (2002) claimed a detec- 
tion of the subordinate line, the derived Li abundances 
are 0.2 — 0.5 dex higher than suggested by the resonance 



line. This is in sharp contrast with our findings. While 
we are not in a position to unequivocally identify the rea- 
sons for the very different conclusions of Ford et al., we 
suspect that the key is the quality of the observations. 
The S/N ratios of the spectra of Ford et al. are between 
200 and 300 around the 610.4 nm line, which may still not 
be sufficient to detect this very weak line in metal-poor 
turn-off stars (note that this line is about 50% stronger 
in a subgiant like HD 140283 with T^s ~ 5700 K than in 
a turn-off star with Tcff « 6200 K for the same Li abun- 
dance). Under these circumstances it is easy to over- 
estimate the line strength due to noise and thus derive 
too high abundances. In comparison, our spectra have 
S/N « 500 or better in the relevant wavelength region. 
Finally, we note that our resolving power is about twice 
as large as obtained by Ford et al. (110000 vs 50 000). 

7. ORIGIN AND EVOLUTION OF LITHIUM 
7.1. Lithium Abundance Trends and Scatter 

Since our sample is relatively small and does not ex- 
tend to extremely low metallicities, our discussion of the 
lithium abundances will be relatively brief. Furthermore, 
our results for the total lithium abundance are in broad 
agreement with many published analyses. The principal 
novelty from our study is an estimate of the isotopic ratio 
for stars on the Spite plateau. 

Our data may be used to quantify the trend of the 
lithium abundance with metallicity (here, Fe/H or 0/H), 
to extrapolate the trend to zero metallicity, and to deter- 
mine the scatter in lithium abundance about the mean 
trend. In applying our data, we have excluded G271- 
162 since no accurate Hct-b ased Teff could be derived 
from this star (see Sec. I4.1f) . Iln addition, the Li- rich 
star HD 106038 has been excluded because it has strong 
overabundances in Si, Ni, Y and Ba (Nissen & Schuster 
1997), and a very high Be abundance relative to other 
halo stars of similar [Fe/H] (Primas et al., in prepara- 
tion). This peculiar abundance pattern suggests that 
the star has undergone a highly unusual nucleosynthesis 
enrichment, possibly due to mass transfer, which justi- 
fies its exclusion here. Such Li-rich stars as HD 106038 
appear to be extremely rare. Perhaps the only known 
other example is the subgiant BD-h23°3912 with [Fe/H] 
= -1.5 and logeu = 2.6 (King et al. 1996). In ad- 
dition, our sample deliberately excluded any of the rare 
Li-deficient halo dwarfs. The nature of these stars is cur- 
rently under debate (cf. Ryan et al. 2002; Frebel et al. 
2005; Aoki et al. 2006; Charbonnel & Primas 2005). Our 
assumption is that these stars have also undergone a very 
special event, a binary merger perhaps (Ryan & Elliott 
2005). 

From the abundances listed in Table 01 the mean 
^Li abundance of the remaining 22 stars is logeT^ — 
2.21 ± 0.07. As shown in Fig. |211 there is, however, a 
pronounced trend with metallicity, which resembles the 
correlation with [Fe/H] found in some previous Li stud- 
ies (e.g. Thorburn 1994; Norris et al. 1994; Ryan et al. 
1996, 1999; Bonifacio et al. 2003, 2006). A univariate 
linear least-square-fit taking into account the estimated 
uncertainties in [Fe/H] and the Li abundances (excluding 
HD 106038 and G 271-162 for the above reasons) implies 

logevLi = (2.409 ±0.020) -I- (0.103 ±0.010) -[Fe/H]. (2) 
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Fig. 22 — Non-LTE "^hi abundances (ID) from the Li I 670.8 nm 
line as a function of metallicity; as for Fig. 1131 solid circles here 
denote stars for which ^Li appears to have been detected while 
open circles are non-detections. In the upper panel logarithmic 
abundances are shown (log CLi vs [Fe /H] ) while in the lower panel 
linear abundances (^Li/H vs Fe/FeQ) are used, which are more ap- 
propriate when extrapolating to the primordial composition. The 
solid lines denote the best linear fits to the data when excluding 
G 271-162 (no Ha-based T^g) and HD 106038 (abnormally high Li 
abundance and other peculiar abundances, see text). In the lower 
panel, the dashed line corresponds to the linear fit when only in- 
cluding the stars with [Fe/H] < -2.2 (G 271-162 again excluded). 
The horizontal solid line in the upper panel shows the predicted 
primordial Li abundance from Big Bang nucleosynthesis and the 
baryon density inferred from WMAP (Spergel et al. 2003) together 
with the corresponding ilcr uncertainties as dotted lines. 



The observed slope is very similar in value to that 
presented by Ryan et al. (1999): 0.103 ± 0.010 vs 
0.118 ± 0.023 for Ryan et al. The scatter around our 
fit is only aohs = 0.033 dex. None of the remaining 22 
stars is identified as an outlier by statistical tests. Re- 
peating the exercise using [0/H] instead of [Fe/H] as the 
independent variable yields 

logevLi = (2.372 ± 0.017) + (0.115 ± 0.011) • [0/H]. (3) 

Here the non-LTE abundances from the O i 777 nm 
triplet have been used to estimate [0/H]. Again, the scat- 
ter is very small: (Tobs — 0.034 dex. 

In Sect. 14.11 we estimated the relative precision in 
our adopted Tcff-values to be about 30 K, which would 
introduce an uncertainty in Li abundance of 0.022 dex. 
Adding also the minor uncertainties from the S/N and 
the other stellar parameters gives a total expected ran- 



dom error of 0.026 dex. This expectation is only slightly 
less than the observed scatter of CTobs — 0.033 dex around 
the linear fit. The indication is that the star-to-star scat- 
ter in lithium abundances at a given [Fe/H] is very small, 
say 0.02 dex or less. (If the temperature error were in- 
creased to 40 K, the expected scatter in the abundances 
would match the observed scatter.) The observed scat- 
ter around the linear fit to our data is very similar to 
that found by Ryan et al., who achieved similarly high 
precision in T^s (ctt ~ 30 K) but utilised lower quality 
spectra {S/N « 100 - 150 and R « 50000). They found 
fobs = 0.031 dex after rejecting one star, BD -t-09°2190, 
as an outlier; we note that with the equivalent width 
measured from our S/N w 650 spectrum, it would no 
longer be considered an outlier. A different and espe- 
cially a larger sample of stars may shed new light on 
the intrinsic scatter in the lithium abundances. For ex- 
ample, Nissen et al. (2005b) find that the otherwise al- 
most identical stars G 64-12 and G 64-37 differ in their Li 
abundances by about 0.16 dex. Their reported observed 
scatter is 0.041 dex when including G 64-12 but 0.035 dex 
when excluding it. No doubt, any intrinsic scatter must 
be very small: < 0.04 dex. 

The near-common belief has been that the primordial 
lithium abundance is obtainable by extrapolation of the 
observed abundances to zero metallicity. One supposes 
that extrapolation should be performed using Fe/H (or 
an alternative such as 0/H) and not the logarithmic 
[Fe/H] (or [0/H]). As seen in Fig. [221 on such a hn- 
ear scale the most metal-poor stars seem to indicate a 
steeper slope than the less metal-poor stars, but there is 
no known explanation for the apparent change of slope 
close to the origin Fe/Fesun = 0. Somewhat different re- 
sults are, therefore, obtained depending on which stars 
are included. When including only the nine stars with 
[Fe/H] < -2.2 (as usual, excluding G 271-162), the esti- 
mated primordial ^Li abundance is ^Li/H « 1.1 ■ 10"^*^ 
or log Eli ~ 2.04. Taking all 22 stars into account results 
in an increase to ^Li/H « 1.4 • 10^^" of logeLi ~ 2.15. 
Very similar results are obtained using extrapolation to 
0/0©= 0. If the presence of ^Li at the level of that 
now observed for LP 815-43 is a general feature among 
the most metal-poor stars and it is produced in spal- 
lation or a-fusion reactions, the ""Li abundances would 
have a minor non-Big Bang contribution which would 
reduce the estimated primordial ^Li abundance further 
by about 0.05 dex (Smith et al. 1998; Ryan et al. 1999; 
2000). 

Our estimates of the primordial abundance do not dif- 
fer significantly from the majority of published determi- 
nations, although they may be somewhat on the low side 
of previous estimates (e.g., Spite & Spite 1982; Thorburn 
1994; Ryan et al. 1996, 1999; Bonifacio & Molaro 1997; 
Bonifacio et al. 2003, 2006). Our estimate of the pri- 
mordial abundance is based on ID model atmospheres 
and includes a correction for non-LTE effects. Presently 
available calculations (Asplund et al. 2003; Barklem et 
al. 2003; also Sec. 15.7(1 suggest that only a minor revision 
will be needed when the abundance analyses are redone 
using 3D models with inclusion of non-LTE effects. 

It has been claimed by Melendez & Ramirez 
(2004) that their new calibration of the infrared flux 
method (IRFM) of effective temperature determinations 
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(Ramirez & Melendez 2005a, b) raises substantially the 
lithium abundances derived from metal-poor dwarf stars. 
This new Teff-scale is significantly hotter at low metal- 
licity than our Ha-based Tcff-scale. It is also hotter 
than the IRFM calibration of Alonso et al. (1996), even 
though the photometric measurements required by the 
IRFM of the two studies agree well for the stars in com- 
mon. For stars with [Fe/H] < -3, such as G 64-12, G64- 
37 and LP 831-70, the new temperatures are 240, 460 
and 320 K, respectively, higher than those estimated us- 
ing the Alonso et al. {V — K) calibration. Ramirez & 
Melendez argue that the strikingly different behaviour 
for the colour calibrations at low [Fe/H] is due to the 
shortage of low metallicity stars in Alonso et al.'s sam- 
ple. 

A detailed scrutiny of Ramirez & Melendez' Toff-scale 
is beyond the scope of the present investigation. Nev- 
ertheless, it is pertinent to estimate the impact of the 
scale on our results. For the most metal-poor stars 
in our sample ([Fe/H] < —2.6), the Tcg-values of Me- 
lendez & Ramirez (2004) are on average 182 ± 72K 
higher than ours based on Ha, which increases the mean 
ID non-LTE ^Li abundance for those particular stars 
from loge7Li= 2.10 to 2.23. Independently of the new 
Tcfi' scale, there are reasons why Melendez & Ramirez 
(2004) arrive at the significantly higher mean ^Li abun- 
dance of log €7 Li = 2.37 ± 0.06 compared with our es- 
timate: their use of convective overshoot model atmo- 
spheres (w -1-0.1 dex), their neglect of non-LTE correc- 
tions (« -1-0.02 dex) and their use of Kurucz models 
rather than MARCS model atmospheres (w -1-0.02 dex). 
The remaining difference with our primordial ''Li abun- 
dance estimate stems from the extrapolation to Fe/Fcgun 
= in the presence of our derived slope with metallic- 
ity. Observational evidence does not favor the convec- 
tive overshoot option as introduced in the Kurucz (1993) 
model atmospheres (e.g. Castelli et al. 1997) nor does 
this approach find any support in 3D hydrodynamical 
simulation of surface convection in metal-poor stars (As- 
plund et al. 1999). Non-LTE abundance corrections, 
while small, should be taken into account. We suggest 
therefore that Melendez & Ramirez' estimate of the pri- 
mordial ^Li abundance should be revised downward by 
at least 0.12 dex regardless of whether or not their Teg- 
scale is preferable to the one employed here. We note 
that our Ha-based temperatures are in excellent agree- 
ment on average with the well-established IRFM {V — K) 
and (& — y) calibrations of Alonso et al. (1996), which 
seems to give little room for such large increases in Tcff as 
advocated by Melendez & Ramirez (2004), In addition, 
accounting for 3D effects should push the Ha-based Teg- 
values downwards. More work is needed to settle the 
Teff-scale of metal-poor stars. 

Recently Bonifacio et al. (2006) have attempted to 
push the abundance determinations to yet lower [Fe/H]. 
In general their results agree well with ours. Their 
mean ID non-LTE Li abundance for 17 stars (exclud- 
ing a subgiant and a spectroscopic binary) is log eli = 
2.17 ± 0.09; they argue that the expected uncertainty 
is about 0.08 dex from the Teff-errors and quality of the 
spectra (median S/N = 180), leaving little or no room 
for intrinsic scatter, in particular when considering the 
possibility of a trend in their Li abundances with metal- 



licity. Due to the small range in [Fe/H] of their sample, 
this correlation may or may not be statistically signifi- 
cant by itself but when including also our results a con- 
sistent picture emerges. We note that their Ha-based 
Teff-scale should be very similar to ours but that their 
[Fe/H] are likely underestimated by 0.1 — 0.2 dex due to 
their use of Fei lines (Sect. 14. 3|) : unfortunately the two 
samples have no stars in common, which prevents a di- 
rect comparison of the parameters and Li abundances. 
Interestingly, the combination of the two samples sug- 
gest a steeper slope for [Fe/H] < —2.3 than for higher 
metallicity. Whether this effect is real or spurious due to 
the analysis methods (ID model atmospheres, Toff-scale 
etc) is an open question for the future. An alternative 
interpretation is that the intrinsic scatter increases sig- 
nificantly for [Fe/H] < -2.3. 

Our major new result concerns the measurements of 
the lithium isotopic ratio and the apparent appearance 
of the analogue for ^Li of the Spite plateau (Fig. IL^ . 
It is impossible to determine if the observed plateau is 
tilted, as the Spite plateau is; as shown in Sect. 17.31 
when accounting for ^Li depletion during the pre-main 
sequence evolution the observed ^Li abundances appar- 
ently increase towards higher metallicity. There may be a 
real scatter associated with this ®Li plateau. The average 
observed abundance from eight stars with a positive ^Li 
detection (HD 106038 is not included) is loge6Li= 0.82, 
i.e. without accounting for any ^Li depletion. Taking the 
mean ^Li/^Li of all stars but HD 106038 (unusually high 
Li abundance) and HD 19445 (expected to have no ^Li) 
together with the mean ^Li abundance of the same stars 
would instead imply log eeLi = 0.54. To both of these esti- 
mates should be added the expected ^Li depletion which 
is > 0.3 dex during the pre-main sequence dur ing which 
^Li remains largely intact, as discussed in Sect. 17.31 If, in 
addition, some ^Li has been destroyed as well this would 
increase the ^Li depletion further. 

7.2. The Lithium Problems 

Since Spite & Spite's (1982) seminal paper, the lithium 
abundance of metal-poor warm dwarfs has been widely 
equated with the abundance provided by the Big Bang. 
The predicted abundance from a standard Big Bang is 
dependent on one free parameter, namely, the baryon 
density. Now, the free parameter is determined inde- 
pendently of any inferred primordial abundances from 
the WMAP-observations of the anisotropics of the cos- 
mic microwave background: ^Li/H — 4.4 ± 1.0 • 10^^° or 
loggYLi = 2.65 ±0.1 (Coc et al. 2004; Cuoco et al. 2004; 
Cyburt 2004). As different calculations yield slightly 
different results even when adopting exactly the same 
baryon density, we have taken an average of the cited re- 
sults. These same calculations show that very little ^Li 
is synthesised, say ^Li/^Li ^ 10^^. 

While the predictions from standard Big Bang nucle- 
osynthesis are in good agreement with the mean deu- 
terium abundance from the six available damped Ly-a 
systems (e.g. Steigman 2005 and references therein), 
they are in substantial conflict with the observations for 
Li. The ^Li problem is now well established: the in- 
ferred primordial ^Li abundance is about 0.5 dex less 
than the prediction. As will become clear from the dis- 
cussion in Sect. 17.41 our survey has in addition identified 
a ^Li problem: the observed high ^Li abundance at very 
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low [Fe/H] can not be explained by neither standard Big 
Bang nucleosynthesis nor apparently by Galactic cosmic 
ray spallation and a- fusion reactions. Possible resolu- 
tions of these problems were reviewed by Lambert (2004). 
It should also be noted that the problem of reconciling 
the inferred and predicted ^Li abundance is not entirely 
decoupled from the goal of achieving an understanding 
of the synthesis of ^Li. Possible modes of ^Li synthesis 
imply simultaneous synthesis of ^Li but this aggravates 
the severity of the ^Li problem. In turn, the presence of 
^Li constrains the possible explanations invoked for the 
^Li problem. 

The solution to the ^Li problem may lie in one or more 
of the following areas: 

Systematic errors in the abundance analysis: The intro- 
duction of an abundance analysis based on 3D model 
atmospheres and non-LTE line formation for the lithium 
resonance and subordinate lines has not resulted in a 
significant upward revision of the lithium abundance ob- 
tained from ID atmospheres and LTE or non-LTE line 
formation. The derived lithium abundance is dependent 
on the adopted Teff-scale but this cannot accommodate 
the higher temperatures necessary to bring inferred and 
predicted abundance into agreement (see Sec. 14.11 and 
I7.1|l . Although the last word has not been said on the 
representation of the physics of the stellar atmospheres 
and the formation of the lithium lines, it would be ex- 
tremely surprising were systematic errors in these areas 
be shown to resolve the lithium problem. 

Present-day abundances do not reflect the primordial Li 
abundance: Dilution and destruction of lithium in metal- 
poor stars are discussed in the following section (Sec. 
[72J). 

Erroneous nuclear reaction rates: Nuclear reactions in- 
cluded in the customary networks have been well stud- 
ied and uncertainties in the reaction rates at the rele- 
vant temperatures for Big Bang nucleosynthesis cannot 
be held responsible for the lithium problem. Coc et al.'s 
(2004) suggestion that addition to the network of a pre- 
viously neglected reaction - ^Be(^D, p)2 "^He - might 
possibly reduce the predicted ^Li abundance is not sub- 
stantiated by a cross-section measurement (Angulo et al. 
2005). (^Be decays to ^Li following Big Bang nucleosyn- 
thesis.) Based on the measured solar neutrino flux, Cy- 
burt et al. (2004) have been able to show that employed 
nuclear data for the ^He(a, 7)^Be reaction can not be in 
error by a large factor as to modify the predicted ^Li 
abundances significantly. 

Non-standard Big Bang nucleosynthesis: Comments on 
the synthesis of the lithium isotopes by non-standard Big 
Bang models are included in our discussion of the syn- 
thesis of ""Li in Sec. 17^ 

7.3. Lithium Destruction and Dilution in Halo Stars 

Identification of the primordial lithium abundance 
from the Spite plateau with the abundance provided by 
the Big Bang rests on the assumption that the lithium 
abundance of the atmosphere of a star on the plateau has 
been unchanged since the star was formed from gas hav- 
ing the lithium abundance provided by the Big Bang. 
Perhaps, the most promising solution to the lithium 
problem represents a failure of this assumption. Solu- 
tions have consequences for the correction of observed 



^Li abundances to the ^Li abundance of the gas from 
which a star formed. 

In standard stellar models, the expected depletion of 
^Li for metal-poor turn-off stars is negligible, i.e., less 
than 0.02 dex (Deliyannis et al. 1990; Pinsonneault et al. 
1992). A larger depletion, > 0.3 dex, is predicted for ^Li 
with most of this occuring during the pre-main sequence 
phase when ^Li is destroyed by protons at the base of 
the star's deep convective envelope (e.g. Richard et al. 
2005). The cross-section for destruction of ^Li is about 60 
times greater than for ''Li. Destruction of ^Li is predicted 
to increase sharply with decreasing stellar mass for main 
sequence stars below the turn-off. Evidently, depletion 
(destruction) of ^Li is according to standard models not 
the answer to the lithium problem but is a factor to be 
considered when we turn to examine synthesis of ^Li. 

The search for a solution to the lithium problem then 
turns to investigation of non-standard models of stellar 
evolution (Lambert 2004). Physical processes neglected 
by standard models but incorporated into non-standard 
models in connection with halo stars in recent years in- 
clude rotationally-induced mixing (e.g. Vauclair 1988; 
Chaboyer & Demarque 1994; Pinsonneault et al. 1999, 
2002; Piau 2005), diffusion^ (e.g. Michaud et al. 1984; 
Salaris & Weiss 2001; Richard et al. 2002; 2005) and 
internal gravity waves (e.g. Talon & Charbonnel 2004; 
Charbonnel & Talon 2005), sometimes working in tan- 
dem with other mechanisms such as mass loss (e.g. Vau- 
clair & Charbonnel 1995) and turbulent mixing. In most 
cases, the predictive power of these processes is restricted 
by their need to introduce one or more free parameters. 

We comment briefly on two non-standard models 
drawn from the suite of proposals: rotationally-induced 
mixing and diffusive separation. An obvious requirement 
of successful non-standard models is that they provide a 
depletion of surface lithium by about 0.5 dex over the 
10-13 Gyr life of the observed stars (we take it as given 
that the stars began life with the WMAP-inferred pre- 
dicted lithium abundance.) Surely, a more demanding 
requirement (and a certain clue to the dominant pro- 
cess) is that this depletion be achieved uniformly over 
the observed sample of stars spanning a range in mass, 
metallicity, age, and rotational angular momentum: the 
depletion sought is about 0.5 dex with a star-to-star scat- 
ter of not more than about 0.03 dex. Our present interest 
is in identifying non-standard models that may solve the 
lithium problem and in obtaining an estimate of the cor- 
rection to be applied to a ^Li abundance to obtain the 
star's original ^Li abundance. 

A description of lithium destruction from rotationally- 
induced mixing depends on processes and initial condi- 
tions that are as yet poorly known, including the initial 
distribution of angular momentum in the star, the trans- 
port of angular momentum through the star and its loss 
from the surface, and attendant mixing of the interior 
material. To make progress requires calibrations using 
observations of rotational velocities and lithium abun- 
dances for stars, e.g., the Sun and stars belonging to 
open clusters. Calibrations have not involved metal-poor 
dwarfs. A probable result of rotational effects on lithium 
depletion is a spread in the final lithium abundance for 

* The term diffusion normally includes the effects of gravitational 
settling, thermal diffusion and radiative acceleration. 
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a sample of stars differing only in their rotational char- 
acteristics. One expects lithium depletion to be most 
severe for the stars initially rotating rapidly and least se- 
vere to vanishingly small for the slowest rotators. Then, 
the observed (small) spread in lithium abundances on the 
Spite plateau is, in principle, an indicator of rotationally- 
induced depletion of lithium. Conversion of the observed 
spread to an indicator of depletion depends on the de- 
scription of rotational-induced mixing and the attendant 
calibrations. 

The recipe and calibrations (Sun, open clusters) in- 
voked by Pinsonneault et al. (1999) show that the de- 
pletion of ^Li {Dj in dex) was related to the dispersion 
(cr in dex) where a/Dy ~ 0.4 and, hence, the required 
depletion by 0.5 implies a = 0.2, a value far in excess of 
the observed upper limit. The depletions of ''Li and ^Li 
were predicted to be correlated: Dj/Dq ~ 0.4 yielding 
k, 1.2 dex. Given that the observed logeeLi ~ 0.8 
for the stars with detected ^Li, the initial abundance 
would have been logeeLj ~ 2.0, the value that has his- 
torically been identified with the Spite plateau for ''Li 
and, more importantly, a value that implies most or all 
of the ''Li came from synthesis by cosmic rays and not 
from the Big Bang! Pinsonneault et al. (2002) argue 
that the Li data of Ryan et al. (1999), who like us mea- 
sured fJobs — 0.03 dex, imply a most-likely ''Li depletion 
of 0.13 dex although destruction factors as high as 0.3 dex 
were not ruled out. The corresponding ^Li destruction is 
then > 0.3 dex, and logeeLj > 1.1 for our stars with ^Li 
detections. In summary, rotationally-induced mixing as 
presented in the models of Pinsonneault et al. appears 
not to be the solution to the ''Li problem. 

Stars on the Spite plateau have a thin convective en- 
velope of uniform composition. Atomic diffusion occurs 
in the radiative zone below the convective envelope. The 
base of the latter mixes with and attains the composi- 
tion of the top of the radiative zone. Diffusion is a slow 
process but is predicted to be effective in the halo dwarfs 
because of their great age and their thin convective en- 
velopes. The diffusive velocity decreases with increasing 
density. Since the stars at the turn-off have the thinnest 
envelopes, the effects of diffusion on surface composition 
are likely to be more severe for these stars than for stars 
further down the main sequence. Lithium is predicted 
to diffuse inwards in the radiative zone and the surface 
abundance to decrease. Inward diffusion of lithium can 
result in its destruction by protons. 

Model calculations including diffusion suggest that the 
lithium abundance for stars on the Spite plateau can 
be reduced to the observed level from a starting value 
corresponding to the WMAP-based prediction (Salaris 
& Weiss 2001; Richard et al. 2002; Richard et al. 
2005). The earhest calculations (Michaud et al. 1984), 
as well as recent ones (Salaris & Weiss 2001), predicted 
lower lithium abundances for stars at the hot end of the 
plateau, in contradiction to the observations. Richard 
and colleagues invoke 'turbulent diffusion' with a recipe 
that restores uniform lithium abundance to the plateau; 
this extra turbulence could for example be induced by 
stellar rotation, which would limit the diffusion but not 
necessarily mix Li to sufficiently high temperatures where 
it is burnt by nuclear reactions as in the models by Pin- 
sonneault et al. (1999, 2002) and others described above. 
Richard et al. (2002, 2005 and private communication) 



have presented detailed calculations of the combined ef- 
fects of diffusion and a parameterized recipe for turbu- 
lent diffusion. The predicted ''Li and ^Li depletion fac- 
tors for our program stars based on their calculations 
are given m Table El We note that the pre-mam se- 
quence depletion factors are likely more uncertain than 
the main sequence and post-main sequence counterparts 
(Proffitt & Michaud 1989). Furthermore, we note that 
the predictions are presently available only for a very 
limited number of metallicities, and thus many depletion 
estimates in Table |H1 are based on extrapolations, which 
likely introduce further uncertainties, in particular to- 
wards higher [Fe/H]. For the disk stars HD 68284 and 
HD 130551, in which ^Li has been detected (Nissen et al. 
1999), we estimate the pre-main sequence ^Li depletion 
to be about 1 dex; these two stars have nearly identi- 
cal mass and [Fe/H], which explains their similarity in 
depletion in spite of their Teg differences. 

The following results are pertinent to our discussion. 
One, there is a recipe for the turbulent diffusion (the so- 
called T6.25 model) that reduces the Big Bang lithium 
abundance by 0.4 — 0.5 dex for turn-off stars to the ob- 
served value and preserves a largely flat plateau (Richard 
et al. 2005). Two, simulations suggest that stars of dif- 
ferent ages and metallicities may have lithium quite uni- 
formly depleted, but a quantitative comparison with the 
small scatter reported by us and others was not made. 
Third, ^Li is partly destroyed during the pre-main se- 
quence evolution even without any accompanying ''Li de- 
struction; due to the lower Toff in general for our higher 
[Fe/H] stars, the observed ^Li plateau becomes tilted 
when accounting for this depletion (Fig. I23|l . Fourth, any 
substantial ''Li destruction is accompanied by an even 
larger ^Li depletion; for example, the T6.25 model which 
provides a 0.4 dex depletion of ''Li predicts 1.6 dex de- 
pletion for ^Li in LP 815-43, in addition to about 0.3 dex 
^Li destruction during the pre-main sequence phase. It 
then follows that the stars with ^Li detections began with 
logevLj = 2.6 and logeeLj « 2.6 too - not a combination 
that can be said to reconcile observations with the Big 
Bang predictions. The choice of a model that provides a 
smaller ''Li depletion also gives a smaller ^Li depletion; 
for example, the T6.09 model reduces ''Li by about 0.25 
dex (lithium problem is lessened but not eliminated) and 
^Li by about 0.6 dex including the pre-main sequence de- 
struction (the initial abundance was logeeLj w 1.4). The 
question is whether alternative recipes for turbulent dif- 
fusion can resolve the ''Li problem without requiring very 
high initial ®Li abundances. 

Although no quantitative predictions yet exist for Pop- 
ulation II stars, a further promising candidate for Li de- 
pletion is internal gravity waves, which lead to angu- 
lar momentum transport and thus mixing (e.g. Garcia 
Lopez & Spruit 1991; Talon & Charbonnel 2004; Char- 
bonnel & Talon 2005). Recent calculations including 
these effects successfully reproduce both the quasi-flat 
solar rotation profile and the degree of Li depletion in 
disk open clusters as a function of age (Charbonnel & 
Talon 2005). We are eagerly awaiting similar modelling 
for metal-poor stars with predictions for the ''Li and ^Li 
depletions, in particular whether large Li depletions can 
be achieved while still retaining the small observed abun- 
dance scatter. 

Our commentary on rotationally-induced mixing and 
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diffusion shows that plausible (rotationally-induced mix- 
ing) and seemingly inevitable (diffusion) processes may 
solve the ^Li problem, but the observers' challenge to 
fit the Spite plateau's shallow slope with respect to T^s 
and [Fe/H] and its smoothness are as yet unmet. These 
potential solutions to the ^Li problem indicate that the 
observed ^Li abundance is almost certainly a lower bound 
to the initial ^Li abundance. Some solutions lead to very 
high initial ^Li abundances that, if correct, imply that 
our understanding of the pre-galactic and Big Bang nu- 
cleosynthesis is poor indeed. 

7.4. Lithium Production in the Early Universe 

The observed ^Li abundance is several orders of mag- 
nitude larger than that predicted from the standard Big 
Bang. On the assumption that the standard Big Bang 
sequence is the correct representation of the primordial 
fireball, synthesis of ^Li is attributed to collisions in the 
intergalactic or interstellar medium between high-energy 
particles (cosmic rays) and ambient nuclei. The two lead- 
ing processes are the a+a fusion reactions and spallation 
reactions involving protons (also, as) and -'^^O (also ■^'^C 
and ^"^N) nuclei. Both processes also produce Li with 
^Li/^Li'- 1 to 2 (Mercer et al. 2001). Only the spalla- 
tion reactions produce ^Be, ^°B, and ^^B. In the classical, 
or direct, Galactic cosmic ray spallation scenario as for- 
mulated by Reeves et al. (1970) and Meneguzzi et al. 
(1971), the accelerated particles are the protons etc hit- 
ting ambient interstellar CNO nuclei. As a consequence, 
this reaction is secondary in nature: [Be/H] oc[0/H] ^. 
In inverse spallation the collisions instead occur through 
fast CNO nuclei and interstellar protons and a-particles. 
Since the CNO nuclei are both produced and acceler- 
ated by supernovae, the resulting spallation yields will 
be essentially independent of the metallicity, a primary 
process: [Be/H] oc[0/H]. One possible way of achieving 
this is through superbubbles (Higdon et al. 1998; Parizot 
& Drury 1999 but see Prantzos 2005, 2006 for a critique 
of this scenario): a sequence of supernovae in rapid con- 
cession in a cluster or OB association accelerates nuclear- 
process enriched material. Establishing the behaviour of 
[0/Fe] with metallicity is clearly important for estimat- 
ing the spallation production rate but unfortunately no 
unanimous verdict have appeared to date in this regard 
(e.g. Israehan et al. 1998, 2001; Boesgaard et al. 1999; 
Asplund & Garcia Perez 2001; Nissen et al. 2002; Ful- 
bright & Johnson 2003; Garcia Perez et al. 2006). 

Since spallation reactions are thought to be the sole 
process by which Be is synthesised and the Be abundance 
is obtainable for stars on the Spite plateau (Boesgaard 
et al. 1999; Primas et al., in preparation), the efficiency 
of a proposed scenario invoking spallation may be rather 
well calibrated using observed Be abundances. The si- 
multaneous production from the a + a fusion reactions 
is not simply calibratable. The observed ^Li/Be ratio 
shows that the ^Li must be produced primarily from the 
fusion reactions at low metallicities with Be necessarily 
from spallation; the meteoritic ^Li/Be ratio (w 6) is a 
representative measure of the production ratios for spal- 
lation but the ratio is 40 for HD 84937 (Smith et al. 
1998) and may reach 150 for LP 815-43 (Primas et 
al., in preparation). The initial ^Li/Bc ratios were likely 
even greater due to ®Li depletion, underscoring the need 
for a + a fusion reactions. 
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Fig. 23. — Upper panel: Observed logarithmic abundances of 
^Li (open triangles) and ^Li (solid circles) as a function of [Fe/H] 
for our program stars; 3a upper limits to the ^Li abundances are 
denoted with arrows. Also shown as open circles are the ^Li detec- 
tions in the halo turn-off star HD 84937 (Smith et al. 1993, 1998; 
Cayrel et al. 1999) and in the two Galactic disk stars HD 68284 
and HD 130551 (Nissen et al. 1999). The large circle correspond to 
the solar system meteoritic ^Li abundance (Asplund et al. 2005), 
while the horizontal solid line is the predicted ^Li abundance from 
Big Bang nucleosynthesis and the baryon density as determined 
by WMAP (Spergel et al. 2003; Coc et al. 2004; Cuoco et al. 
2004; Cyburt 2004). The (lower) solid, dashed, (lower) dotted and 
dash-dotted lines correspond to the models for cosmic ray ''Li pro- 
duction by Prantzos (2006), Ramaty et al. (2000), Fields & Olive 
(1999) and Vangioni-Flam et al. (2000), respectively; see text for 
details. The higher dotted line shows the ^Li abundance assum- 
ing the Fields & Olive model with a ^Li/^Li production ratio of 
1:1.5 and a primordial abundance of loge7Li= 2.1. Lower panel: 
Same as in the upper panel but taking into account the predicted 
pre-main sequence ®Li depletion from Richard et al. (2002, 2005; 
and 2005, private communication) to the observed ''Li abundances; 
these are the minimum values corresponding to negligible ^Li de- 
pletion (see Sect. 17. 31 for details). 



A large number of Galactic chemical evolution mod- 
els including spallation and a-fusion reactions have been 
presented in recent years and compared with observed Li, 
Be and B isotopic abundances (e.g. Yoshii et al. 1997; 
Fields & Olive 1999; Parizot 2000; Vangioni-Flam et al. 
2000; Ramaty et al. 2000; Suzuki & Yoshii 2001; Al- 
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ibes et al. 2002; Prantzos 2005, 2006). Four representa- 
tive predictions are presented in Fig. ESI the model by 
Fields & Olive (1999) is based on direct cosmic ray spal- 
lation but with a very steep [0/Fe] trend (Israelian et 
al. 1998; Boesgaard et al. 1999); the model by Ramaty 
et al. (2000) is for inverse spallation with special con- 
sideration of the energetics of the light element produc- 
tion; the Vangioni-Flam et al. (2000) model is for inverse 
spallation in supcrbubbles with an assumed low energy 
(< 30 MeV/nucleon) cosmic ray component (Meneguzzi 
& Reeves 1975); and the calculations by Prantzos (2006) 
assume that the cosmic ray composition does not vary 
with time and that the light isotopes are produced as pri- 
maries, i.e. like in the inverse spallation models. These 
calculations fit about equally well the observed run of the 
Be abundance with the Fe abundance, although Ramaty 
et al. (2000) argue that the Fields & Olive Be produc- 
tion is significantly overestimated through the adoption 
of unrealistic Fe yields . A fair fit to the observed ^Li 
abundance in HD 84937 (previously the sole data point 
for this isotope) when ignoring the inevitable pre-main 
sequence destruction is achieved for the Fields & Olive 
and Vangioni-Flam et al. models but it requires some as- 
sumptions to maximise ^Li production. There are other 
problems and difficulties with the calculations: for ex- 
ample, the efficiency assumed for cosmic ray acceleration 
may be unrealistic given the available kinetic energy from 
supernovae, which may rule out the low energy cosmic 
ray scenario (Ramaty et al. 2000; Suzuki & Yoshii 2001; 
Prantzos 2005, 2006). Furthermore, several of the spalla- 
tion models presented here invoke a [0/Fe] trend that is 
significantly higher than observed by Nissen et al. (2002) 
and us (Sect. 14.4(1 : note that in addition our adopted so- 
lar O abundance is 0.1 — 0.2 dex lower than assumed by 
Fields & Olive and Ramaty et al. As a consequence, their 
predicted ^Li abundances should be shifted downwards 
by « 0.3 dex for [Fe/H] < —2. Finally, the subsequent 
downward revision by a factor of two of the measured 
cross-section for ^Li production by a-fusion (Mercer et 
al. 2001) will further push the predictions to lower yields 
for all but the Prantzos model. 

Now that additional ^Li detections are available with 
the high ^Li abundance in LP 815-43 in particular, it 
is clear that the model calculations do not reproduce 
the ^Li measurements and the apparent observed ^Li 
plateau. Of the here presented models, the one by Prant- 
zos (2006, shown as a solid line in Fig. I23|l is likely 
the most realistic as it requires the least number of as- 
sumptions in terms of maximising the production effi- 
ciency and is based on the most recent reaction cross- 
sections. Still we note that it makes the as yet unproven 
assumption that the cosmic ray composition stays con- 
stant with time and thus may well overestimate the ^Li 
production. The discrepancy between models and obser- 
vations at low [Fe/H] is further aggravated when account- 
ing for the pre-main sequence ^Li depletion of > 0.3 dex 
(Table EJ, as seen in Fig. invoking any significant 
amount of ""Li depletion would make the situation worse 
still (Sect. I7.3fl . The inclusion of ^Li pre-main sequence 
depletion tilts the ^Li observations to follow roughly the 
metallicity-dependence of the models, albeit at a higher 
^Li abundance; we recall here that the pre-main sequence 
depletion calculations are rather uncertain and that most 
of the ^Li depletion may in fact take place during the 



main sequence. At face value, it seems like a remark- 
able coincidence that the ^Li production from cosmic ray 
spallation and Li-depletion can result in a nearly flat ^Li 
abundance plateau as observed. However, we do note 
that we are only able to detect ^Li/^Li ratios of > 0.02 
for observational reasons and thus would not be able to 
identify the stars in which the depletion is much larger 
than the corresponding production. 

A means to boost the light element production over di- 
rect and indirect cosmic ray spallation in the interstellar 
medium has recently been put forward by Nakamura et 
al. (2006, see also Fields et al. 1996 and Nakamura & 
Shigeyama 2004) . They speculate that the fast CNO-rich 
ejecta of very energetic (« 10^^ erg — 10^^ J) supernovae 
having lost all but a small fraction of their H- and He- 
rich envelopes can interact directly with the circumstellar 
material without shock acceleration of cosmic rays. Such 
SNe Ic events may produce significant amounts of pri- 
mary ^Li and ^Be under suitable conditions. Nakamura 
et al. (2006) have identified a parameter space that seem 
to explain the high ^Li/0 ratio we observe in LP 815-43 
but it requires some fine-tuning in regards the explosion 
energy and the amount of available He in the ejecta. It 
remains to be seen whether real supernovae fulfill these 
requirements. 

An additional and alternative scenario places lithium 
synthesis by a -t- a fusion reactions during the forma- 
tion of the Galaxy within the framework of hierarchical 
structure formation (Suzuki & Inoue 2002, 2004). In- 
falling and merging subgalactic clumps would have pro- 
duced shocks that could have accelerated the a-particles 
at early epochs. Due to our incomplete understanding 
of the details of galaxy formation, models for the ac- 
companying ^Li production by necessity include several 
free parameters. To achieve the observed high ^Li abun- 
dances at [Fe/H] < — 2.3 requires strong structure forma- 
tion shocks to have been in place very early on (k, 0.1 Gyr 
after Big Bang) but vanished relatively shortly there- 
after (« 0.1 Gyr). An argument against this scenario 
is that mergers in ACDM continue to the present day 
with only a relatively small dark matter halo present at 
early times, which would appear to make it difficult for 
the necessary strong shocks to have formed during the 
required time interval (Prantzos 2006). RoUinde et al 
(2005, see also Montmerle 1977) have presented a simi- 
lar scenario of pre-Galactic production of ^Li (and ^Li) 
due to an ad-hoc early burst of cosmological cosmic rays, 
possibly due to Population HI stars. The energetics of 
such a production were not investigated. Prantzos (2006) 
also speculate that accretion to a supermassive black hole 
such as the one believed to exist currently in the Galactic 
center could generate sufficient energy to accelerate the 
necessary cosmic rays. This scenario requires some ex- 
treme assumptions about the time of its formation and 
its efficiency in converting accretion energy to particle 
acceleration to produce sufficient ^Li at [Fe/H] « —2.7. 

The final potential site for ^Li synthesis is the Big 
Bang. In a standard Big Bang much too little ^Li is 
produced to explain our observations. An extension of 
the standard model for particle physics such as super- 
symmetry predicts the existence of various exotic parti- 
cles, including the gravitino, axion and the neutralino; 
the first two of these particles are predicted to be un- 
stable but the neutralino is likely stable. The decay and 
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the annihilation during or shortly after the era of Big 
Bang nucleosynthesis can alter the resulting light ele- 
ment abundances, provided the masses and life-times of 
these putative particles are right (e.g. Dimopoulos et 
al. 1988; Jedamzik 2000, 2004a,b, 2005; Kawasaki et al. 
2004; Ichikawa et al. 2004). The annihilation of the neu- 
tralino can release sufficient energy to produce ^Li by 
non-thermal reactions like '^H(Q;,n)^Li and '^He(a,p)^Li 
(Jedamzik 2004a). The amount of ^Li produced depends 
on the mass of the neutralino and the exact annihilation 
channels but can reach the levels observed in some of 
our very metal-poor stars; this scenario does not, how- 
ever, resolve the ^Li dilemma described above. Another 
option to produce ^Li is through the decay of particles 
like the gravitino and axion. While the electromagnetic 
decay routes of such particles result in ■^He/D ratios in- 
consistent with observations (Ellis et al. 2005), the in- 
jection of energetic nucleons through the hadronic decay 
about 10^ s after Big Bang can lead to substantial ^Li 
production without spoiling the agreement with D and 
the He isotopes. Indeed, for the right combination of par- 
ticle properties, a simultaneous production of ^Li and de- 
struction of ^Li appears achievable and may explain both 
the observed ^Li plateau and the low ^Li abundances 
in comparison with standard Big Bang nucleosynthesis 
(Jedamzik 2004b, 2005). Thus, both of the Li problems 
can conceivably be solved at the same time. While these 
idea are very attractive, they rests on as yet unproven 
and speculative physics. 

8. CONCLUDING REMARKS 

Remarkable as it may have seemed in 1982 following 
the discovery of the Spite plateau and the near-universal 
association of the plateau's lithium abundance with pri- 
mordial nucleosynthesis, the debate over that association 
has not been stilled. Two key recent factors heating the 
debate have been the measurement of the Big Bang's 
photon-to-baryon ratio (77) using the WMAP measure- 
ments of the anisotropics in the cosmic microwave back- 
ground and the detections, here and previously, of the 
lighter isotope ^Li. A standard model of the Big Bang 
with the WMAP estimate of predicts a lithium (^Li) 
abundance greater by about 0.5 dex than the value gen- 
erally attributed to the Spite plateau, and a ^Li abun- 
dance several orders of magnitude below the value re- 
ported here. 

These two discrepancies resist simple resolution, as 
should be clear from the discussion in Section 7. One sup- 
poses that resolution of the discrepancies may be found 
in one or more of the following three areas: 
1. Composition of the gas from which stars on the Spite 
plateau formed 

Ignoring the possibility that the ^Li was synthesised af- 
ter the formation of the star, a conjecture is that ^Li 
was synthesised either in a non-standard Big Bang (e.g., 
hadronic decays of exotic particles resulting in post- 
processing the ^Li from a standard Big Bang) or through 
a — a fusion reactions in the proto-Galaxy (e.g. hierar- 
chical structure formation shocks and/or supermassive 
black holes). Introduction of a non-standard Big Bang 
may lower the predicted ^Li abundance and increase the 
predicted ^Li abundance. Qualitatively, these changes 
resolve both discrepancies provided that points 2 and 3 
below do not come into play. Addition of lithium synthe- 



sis from fusion reactions in pre-Galactic g the sole 
additional site over a standard Big Bang, may remove 
the ®Li discrepancy but may aggravate the discrepancy 
between the predicted ^Li abundance and the observed 
abundance, unless one adopts the heresy that the WMAP 
estimate of rj should be replaced by the lower value that 
results in the minimum synthesis of ''Li. (The nuclear 
reaction network appears to be free of significant uncer- 
tainties.) 

2. Relationship between the stars ' initial composition set 
10 - 13 Gyr ago and the present composition of their 
stellar atmospheres 

There is tension on this point between those who stress 
the significance of the observed very small scatter of the 
lithium abundances along the Spite plateau and those 
who identify physical processes that ought to change 
the atmospheric compositions of plateau residents (Sec- 
tion 7.4). The very small scatter would seem to im- 
ply that the change in surface lithium abundance has 
been itself slight. Yet, some argue that reduction of 
the lithium abundance through diffusion, perhaps sup- 
plemented by rotationally-induced mixing and other pro- 
cesses, must occur but can result in a small star-to-star 
scatter of lithium abundances. The presence of ^Li, as 
we noted, complicates interpretations that rely on de- 
pletion/destruction of "^Li to reduce the WMAP-based 
Big Bang ^Li abundance to the observed value for the 
Spite plateau; very high to extraordinarily high initial 
^Li abundances may be implied. There are clearly open 
questions for pursuit by observers. 

Diffusion does not exclusively act on lithium. Rela- 
tive abundances of abundant elements are changed but 
these elements are not subject to destruction as lithium 
is. Alterations of the abundances are lessened and then 
removed as stars evolve from the main sequence to sub- 
giants and then to giants. Careful comparisons of stars 
along this sequence might confirm that diffusion was op- 
erable in their main sequence progenitors. At present, 
abundance data from dwarfs and giants have given seem- 
ingly consistent results for composition as a function of 
metallicity, i.e., the same run of [X/Fe] versus [Fe/H] 
(e.g. Gratton et al. 2001). 

3. Measurement of the present ^Li and Li abundances 
Systematic errors are now unlikely to be uncovered that 
will lift the Spite plateau's lithium abundance by the 0.5 
dex to match the Big Bang WMAP-based abundance 
(Section 7.1). The adjustments (e.g, the necessary in- 
crease of effective temperatures by about 700 K) unless 
specific to lithium would result in large changes to abun- 
dances of other elements. Presently, there is generally 
good agreement between the compositions of metal-poor 
dwarfs and giants of the same metallicity and this agree- 
ment would be compromised unless adjustments were 
also made to the abundance analyses of giants with their 
very different atmospheric structures. The simplest con- 
ceivable adjustment - raising the effective temperatures 
by 700 K - introduces more problems than the one it 
solves. 

Absence of systematic errors affecting the determina- 
ton of the ^Li/^Li ratio cannot be so strongly dismissed: 
the presence of ^Li is deduced from a slight additional 
asymmetry to the intrinsically asymmetric ''Li line. The 
various tests outlined in Sect. |S1 give confidence that 
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^Li has indeed been detected. This confidence should be 
tested by observations of additional stars like HD 19445 
for which pre-main sequence destruction of ^Li is so great 
that this isotope should be undetectable. Obviously, ad- 
ditional ^Li observations at very low [Fe/H] would be 
extremely valuable, including a confirmation of our de- 
tection in LP 815-43. 

The confidence commonly and frequently expressed af- 
ter 1982 that primordial lithium had been discovered and 
measured and could be matched with a prediction from a 
standard Big Bang is today weakened. There remains a 
belief that much of the lithium is traceable to primordial 
nucleosynthesis but quantitative agreement between pre- 
diction and observation of the ®Li and ^Li abundances in 
stars on the Spite plateau is proving elusive. Observers 
will suppose that well designed additional spectroscopic 
analyses of metal-poor stars will provide some of the clues 
necessary to determining the primordial abundance of 
the two lithium isotopes. Theoreticians may combine 
these clues with other evidence and probe deeper into 
primordial nucleosynthesis. Interest in the lithium of the 



Spite plateau is greater than in 1982, even as the lack 
of concordance between predictions of the standard Big 
Bang and measured lithium abundances has been demon- 
strated. 
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Table 1. Observational details, including date of observation, Univer- 
sal Time for the middle of the exposure, total exposure time, signal- 
to-noise AROUND the Li I 670.8 nm line in the reduced spectrum, and the 
measured heliocentric radial velocity of the star. 



otar 


V 


Date 


UTC 


^exp 


S/N 


frad 




[mag] 






[min] 




n — 1 1 

[kms ^\ 


11 LJ OOU 1 


9.26 


9000 07 94 


08:43 


45 


560 


—47.3 


ILiy _L r?' I'lfJ 


8.05 


2000-08-04 


10:13 


20 


790 


— 139.7 


HD 59392 


9.73 


2002-02-05 


02:53 


60 


710 


267.9 


HD 1022(10 


8.74 


2002-02-05 


07:36 


40 


610 


161.1 


HD lOfiOSR 

L1.1-/ A-\J\j\JkJ\J 


10.18 


2000-07-23 


23:57 


60 


450 


99.6 






2000-07-24 


23:56 


40 


350 


99.9 




7.21 


2000-07-24 


23:15 


13 


1020 


— 170.4 






2002-02-05 


08:52 


15 


770 


— 170.4 






2on'i-ns-2s 


00:38 


19 


820 


— 170.6 


HD IfiOfil? 


8.73 


2000-07-25 


01:34 


25 


650 


100.6 




9.66 


2000-08-04 


07:23 


45 


640 


50.8 




10.05 


2002-02-05 


06:24 


80 


470 


198.3 


HD 338529 


9.37 


2000-07-25 


04:12 


60 


680 


— 128.3 




10.00 


9000-07-9'=; 

^UUU U 1 


00:43 


60 


530 


58.4 


n 020-024 


11.13 


2000-07-24 


02:19 


90 


430 


34.2 






9000-07-9^ 


02:47 


90 


420 


34.5 


n 07^-n'^i 

vjr \J 1 'J yjiji. 


10.52 


2000-07-24 


09:46 


60 


500 


57.8 






2000-08-04 


09:26 


60 


500 


57.8 




9.48 


9000 0^^ 04 


06:30 


45 


700 


—288.6 




11.15 


2002-02-05 


04:25 


90 


400 


266.6 






2002-02-06 


05:16 


90 


400 


266.0 


BD +03° 0740 


9.80 


2002-02-05 


01:03 


60 


740 


174.2 


BD- 13° 3442 


10.29 


2002-02-06 


06:27 


90 


400 


116.1 


CD -30° 18140 


9.95 


2000-07-24 


05:51 


60 


450 


17.9 






2000-07-25 


06:03 


25 


350 


17.9 


CD -33° 1173 


10.94 


2002-02-06 


01:13 


90 


610 


47.6 


CD -33° 3337 


9.08 


2002-02-05 


01:57 


30 


740 


73.7 


CD -35° 14849 


10.57 


2000-07-24 


07:27 


90 


550 


108.0 






2000-07-25 


06:33 


25 


350 


108.1 


CD -48° 2445 


10.55 


2002-02-06 


03:01 


100 


640 


319.2 


LP 815-43 


10.91 


2000-07-24 


04:12 


90 


400 


-3.8 






2000-07-25 


05:35 


40 


300 


-3.5 






2004^08-30 


04:52 


150 


420 


-3.6 






2004^08-31 


01:35 


150 


490 


-3.9 
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Table 2. Derived values of effective temperature, surface gravity, 

METALLICITY (AS MEASURED BY Fe II LINES), OXYGEN ABUNDANCE (NON-LTE 
VALUES FROM THE Ol TRIPLET), MICROTURBULENCE AND ABSOLUTE MAGNI- 
TUDE. 



btar 




logs 


[Fe/HJ 


[O/HJ 


?turb 


My 




[K] 


[cgs] 






n —11 






6026 


4.08 


— 1.14 


—0.79 


1.5 


4.01 




5980 


4.42 


—2.02 


— 1.38 


1.3 


5.13 


HD 59392 


5936 


3.99 


— 1.61 


— 1.04 


1.5 


3.97 


ffn 1 02200 


6062 


4.15 


— 1.25 


—0.85 


1.3 


4.27 


HD 1 OfiO'HS 


5905 


4.30 


— 1.35 


—0.76 


1.2 


4.84 




5753 


3.70 


—2.40 


— 1.74 


1.5 


3.35 


HD 1 fiOfil 7 


5990 


3.79 


— 1.76 


— 1.43 


1.5 


3.35 


HD 213657 


6180 


3.92 


— 1.90 


— 1.34 


1.5 


3.60 


HD 29S9Sfi 


6103 


4.22 


— 1.33 


—0.98 


1.3 


4.41 


HD 33S529 


6335 


4.04 


—2.26 


— 1.68 


1.5 


3.73 




6298 


3.99 


—2.30 


—1.71 


1.5 


3.70 


G 020-024 


6247 


3.98 


-1.89 


-1.35 


1.5 


3.72 


G 075-031 


6000 


4.08 


-1.02 


-0.68 


1.3 


3.96 


G 126-062 


6183 


4.11 


-1.51 


-1.05 


1.5 


4.05 


G 271-162 


6230 


3.93 


-2.30 


-1.79 


1.5 


3.61 


BD -|-09°2190 


6392 


4.09 


-2.66 


-2.22 


1.5 


3.89 


BD +03°0740 


6266 


4.04 


-2.65 


-2.06 


1.5 


3.93 


BD -13°3442 


6311 


3.86 


-2.71 


-2.09 


1.5 


3.43 


CD -30° 18140 


6222 


4.11 


-1.90 


-1.28 


1.5 


4.06 


CD -33° 1173 


6390 


4.28 


-2.92 


-2.33 


1.5 


4.36 


CD-33°3337 


5897 


4.01 


-1.31 


-0.76 


1.3 


3.92 


CD -35° 14849 


6244 


4.30 


-2.27 


-1.70 


1.5 


4.56 


CD -48° 2445 


6222 


4.25 


-1.93 


-1.40 


1.5 


4.42 


LP 815-43 


6400 


4.17 


-2.74 


-2.20 


1.5 


4.10 
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Table 3. Li isotopic abundances derived from the Lii 670.8 nm reso- 
nance LINE TOGETHER WITH THE Ll ABUNDANCES ESTIMATED FROM THE Ll I 
610.4 NM SUBORDINATE LINE. BOTH THE ID LTE AND NON-LTE Ll ABUN- 
DANCES ARE LISTED FOR THE TWO LINES. THE QUOTED UNCERTAINTIES IN 
THE LTE ABUNDANCES ARE THOSE ONLY DUE TO THE FINITE S/N AND DO 
NOT FOR EXAMPLE TAKE INTO ACCOUNT ERRORS IN THE STELLAR PARAMETERS. 
The SECOND AND THIRD COLUMN GIVE THE FRO-IECTED ROTATIONAL VELOC- 
ITY -Urotsini (ASSUMED TO BE 0.5 KM S^^ FOR ALL STARS) AND THE MACRO- 
TURBULENCE Cmacro; NOTE THAT INSTRUMENTAL BROADENINC; HAS ALSO BEEN 
APPLIED. The LISTED EQUIVALENT WIDTHS FOR, THE 670.8 ,'\ND 610.4 NM LINES 
ARE THOSE COMING FROM THE BEST-FITTING THEORETICAL LINE PROFILES TO 
THE OBSERVED PROFILES. 



Li I 670.8 nm Li i 610.4 nm 



Star Vrot sin i Cmacro Wx log ETLi log ETLi ^Li/'^'Li \\ \ log ELi log eLi 

[kiris-1] [kiiis-1] [pill] (LTE) (iioii-LTE) [])iii] (LTE) (uou-LTE) 



HD 3567 


0. 


,50 


4, 


.71 


± 


0, 


,19 


3, 


,83 


2, 


.32 


± 


0, 


,01 


2, 


.32 


0.017 


± 


0.012 


0.19 


2, 


.30 


± 


0.04 


2.36 


HD 19445 


0. 


,50 


3, 


.83 


± 


0. 


,13 


3. 


,51 


2, 


.21 


± 


0. 


,01 


2, 


.20 


0.002 


± 


0.008 


0.16 


2, 


.19 


± 


0.05 


2.24 


HD 59392 


0. 


50 


4, 


.37 


± 


0. 


.35 


3. 


.92 


2, 


.24 


± 


0. 


.01 


2, 


.24 


0.021 


± 


0.016 


0.15 


2, 


.16 


± 


0.03 


2.22 


HD 102200 


0. 


50 


4, 


.44 


± 


0. 


.24 


3. 


.30 


2, 


.25 


± 


0. 


.01 


2, 


.25 


0.047 


± 


0.013 


0.18 


2, 


.30 


± 


0.03 


2.36 


HD 106038 


0. 


50 


3 


.88 


± 


0. 


,12 


6. 


.31 


2, 


.49 


± 


0. 


.01 


2, 


.48 


0.031 


± 


0.006 


0.30 


2, 


.45 


± 


0.02 


2.52 


HD 140283 


0. 


50 


4, 


.87 


± 


0. 


.07 


4. 


,78 


2, 


.17 


± 


0. 


.01 


2, 


.20 


0.008 


± 


0.006 


0.18 


2, 


.11 


± 


0.03 


2.20 


HD 160617 


0. 


50 


4, 


.54 


± 


0. 


.18 


4. 


,05 


2, 


.28 


± 


0. 


,01 


2, 


.28 


0.036 


± 


0.010 


0.19 


2, 


.28 


± 


0.03 


2.34 


HD 213657 


0. 


,50 


5, 


.02 


± 


0. 


,08 


3. 


,00 


2, 


.27 




0. 


,01 


2 


.25 


0.011 


± 


0.011 


0.15 


2, 


.24 


± 


0.04 


2.29 


HD 298986 


0. 


50 


4, 


.70 


± 


0, 


,20 


3, 


,04 


2, 


.25 


± 


0, 


,01 


2, 


.24 


0.017 


± 


0.015 


0.15 


2, 


.24 


± 


0.04 


2.29 


HD 338529 


0. 


50 


5, 


.65 


± 


0, 


,09 


2, 


,37 


2, 


.25 


± 


0, 


,01 


2, 


.22 


0.010 


± 


0.013 


0.13 


2, 


.24 




0.05 


2.28 


G 013-009 


0. 


50 


5, 


.85 


± 


0, 


,07 


2, 


,39 


2, 


.22 


± 


0, 


,01 


2, 


.19 


0.048 


± 


0.019 


0.11 


2, 


.15 


± 


0.06 


2.20 


G 020-024 


0. 


50 


5, 


.12 


± 


0, 


,22 


2, 


,49 


2, 


.21 


± 


0, 


,01 


2, 


.18 


0.070 


± 


0.017 


0.11 


2, 


.14 


± 


0.07 


2.18 


G 075-031 


0. 


50 


4, 


.95 


± 


0, 


,22 


3, 


,78 


2, 


.30 


± 


0, 


,01 


2, 


.30 


0.015 




0.012 


0.24 


2, 


.40 


± 


0.03 


2.46 


G 126-062 


0. 


50 


4, 


.76 


± 


0, 


,30 


2, 


,76 


2, 


.26 


± 


0, 


,01 


2, 


.24 


0.001 


± 


0.015 


0.14 


2, 


.22 


± 


0.04 


2.27 


G 271-162 


0. 


50 


5, 


.78 


± 


0, 


13 


2, 


,91 


2, 


.27 


± 


0, 


,01 


2, 


.25 


0.019 


± 


0.012 


0.11 


2, 


.14 


± 


0.05 


2.19 


BD+09°2190 


0. 


50 


6, 


.27 


± 


0, 


,35 


1, 


,66 


2, 


.13 


± 


0, 


,01 


2, 


.10 


-0.033 ±0.027 


< 0.11 




< 


2. 


20 


< 2.25 


BD+03°0740 


0. 


50 


6, 


.71 


± 


0, 


15 


2, 


,08 


2, 


.12 




0, 


,01 


2, 


.11 


0.015 


± 


0.017 


< 0.10 




< 


2. 


10 


< 2.15 


BD-13°3442 


0. 


50 


6 


.38 


± 


0, 


,20 


2, 


,04 


2, 


.16 


± 


0, 


,01 


2 


.14 


0.001 


± 


0.028 


0.11 


2 


.15 


± 


0.08 


2.19 


GD -30° 18140 


0. 


,50 


5, 


.07 


± 


0, 


,11 


2, 


,76 


2, 


.25 


± 


0, 


,01 


2, 


.22 


0.042 


± 


0.013 


0.14 


2, 


.23 


± 


0.03 


2.27 


CD -33° 1173 


0. 


50 


5, 


.82 


± 


0. 


.61 


1. 


.60 


2, 


.08 


± 


0. 


.01 


2, 


.05 


0.013 


± 


0.041 


0.07 


2, 


.01 


± 


0.10 


2.04 


CD -33°3337 


0. 


50 


4, 


.91 


± 


0. 


.14 


3. 


.90 


2, 


.23 


± 


0. 


.01 


2, 


.24 


0.020 


± 


0.010 


0.20 


2, 


.27 


± 


0.03 


2.34 


CD -35° 14849 


0. 


50 


4, 


.68 


± 


0. 


.12 


2. 


,84 


2, 


.27 


± 


0. 


,01 


2, 


.24 


0.020 


± 


0.012 


0.15 


2, 


.26 


± 


0.05 


2.30 


CD -48° 2445 


0. 


,50 


5, 


.18 


± 


0, 


,10 


2, 


,56 


2, 


.22 


± 


0, 


,01 


2, 


.19 


0.032 


± 


0.014 


0.15 


2, 


.27 


± 


0.03 


2.31 


LP 815-4.'^ 


0. 


50 


6 


,20 


± 


0, 


19 


1, 


,89 


2, 


,16 


± 


0, 


01 


2, 


,13 


0.016 


± 


0.022 


0.11 


2, 


,17 


± 


0.06 


2.21 
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Table 4. Effects of the derived ^Li/'^'Li ratio to the stellar parame- 
ters AND SUITE OF ID MODEL ATMOSPHERES 



Star ^Li/'^Li 





norma,! 




Alogp 


nj 


^?turb 


i^urucz* 


HD3567 


0.017 


0.016 


0.015 


0.017 


0.015 


0.018 


HD 19445 


0.002 


0.000 


0.002 


0.002 


0.001 


0.003 


HD59392 


0.021 


0.020 


0.021 


0.021 


0.020 


0.021 


HD102200 


0.047 


0.043 


0.045 


0.047 


0.046 


0.044 


HD106038 


0.031 


0.028 


0.030 


0.032 


0.029 


0.028 


HD140283 


0.008 


0.008 


0.008 


0.008 


0.007 


0.006 


HD160617 


0.036 


0.033 


0.034 


0.035 


0.035 


0.035 


HD213657 


0.011 


0.008 


0.011 


0.011 


0.010 


0.011 


HD298986 


0.017 


0.014 


0.016 


0.017 


0.017 


0.018 


HD338529 


0.010 


0.008 


0.011 


0.010 


0.011 


0.011 


G013-009 


0.048 


0.047 


0.049 


0.048 


0.049 


0.049 


G020-024 


0.070 


0.070 


0.069 


0.070 


0.069 


0.070 


G075-031 


0.015 


0.013 


0.014 


0.015 


0.013 


0.016 


G126-062 


0.001 


-0.001 


0.000 


0.000 


-0.001 


0.003 


G271-162 


0.019 


0.018 


0.019 


0.018 


0.018 


0.019 


BD -|-09°2190 


-0.033 


-0.035 


-0.033 


-0.033 


-0.032 


-0.033 


BD +03° 0740 


0.015 


0.013 


0.015 


0.014 


0.013 


0.015 


BD- 13° 3442 


0.001 


-0.001 


0.001 


0.001 


0.001 


0.001 


CD -30° 18140 


0.042 


0.039 


0.043 


0.041 


0.041 


0.044 


CD -33° 1173 


0.013 


0.010 


0.012 


0.013 


0.013 


0.011 


CD -33° 3337 


0.020 


0.019 


0.019 


0.021 


0.017 


0.021 


CD -35° 14849 


0.020 


0.017 


0.020 


0.020 


0.020 


0.020 


CD -48° 2445 


0.032 


0.031 


0.033 


0.032 


0.033 


0.034 


LP815-43 


0.046 


0.045 


0.046 


0.048 


0.046 


0.046 



''Derived ®Li/^Li when using Kurucz model atmospheres (without convective overshoot as computed by CasteUi et al. 1997) instead of 
the standard MARCS models otherwise used 



Table 5. Estimated projected stellar rotational velocities and 

^Ll/'^'Ll BASED ON A 3D ANALYSIS 



Star 


Viot sin i 


6Li/'^Li(3D) 


6Li/'^Li(lD) 




[kms^-*-] 







HD3567 


0, 


,35 


± 





.60 


0, 


,000 


dz 


0, 


,012 


0, 


,017 




0, 


,012 


HD19445 


2, 


,55 


± 


0, 


,14 


0, 


,030 


± 


0, 


,008 


0, 


,002 




0, 


,008 


HD140283 


2, 


,01 


± 


0, 


,14 


0, 


,023 




0, 


,006 


0, 


,008 


± 


0, 


,006 


HD160617 


1, 


,19 


± 


0, 


,61 


0, 


,052 




0, 


,012 


0, 


,036 


± 


0, 


,010 


HD213657 


1, 


,93 


± 


0, 


.23 


0, 


,030 


± 


0, 


,012 


0, 


Oil 


± 


0, 


,011 


HD338529 


2. 


,50 


± 


0, 


.18 


0, 


,017 


± 


0, 


,013 


0. 


,010 


± 


0, 


,013 


G013-009 


2. 


,77 


± 


0, 


.24 


0. 


,056 


± 


0, 


,021 


0. 


,048 


± 


0, 


,019 


G020-024 


2. 


,30 


± 


0, 


.17 


0. 


,081 


± 


0, 


,016 


0. 


,070 


± 


0, 


,017 


G075-031 


0. 


,67 


± 


0, 


.82 


0. 


,000 


± 


0, 


,012 


0. 


,015 


± 


0, 


,012 


G271-162 


2. 


,65 


± 


0, 


.20 


0. 


,026 


± 


0, 


,012 


0. 


,019 


± 


0, 


,012 


BD +09°2190 


2. 


,84 


± 





.58 


0, 


,000 


± 


0, 


,027 


-0, 


,033 


± 


0, 


,027 


BD+03°0740 


3, 


,80 


± 


0, 


,27 


0, 


,026 


± 


0, 


,019 


0, 


,015 


± 


0, 


,017 


BD-13°3442 


3, 


,02 




0, 


,39 


0, 


,000 


± 


0, 


,028 


0, 


,001 


± 


0, 


,028 


CD-30°18140 


2, 


,21 


± 


0, 


,26 


0, 


,058 


± 


0, 


,014 


0, 


,042 


± 


0, 


,013 


CD-33°1173 


1, 


,81 


± 


1, 


,45 


0, 


,012 


± 


0, 


,027 


0, 


,013 


± 


0, 


,041 


CD -35° 14849 


0, 


10 


± 


0, 


,19 


0, 


,026 




0, 


,011 


0, 


,020 


± 


0, 


,012 


CD -48° 2445 


2, 


,34 


± 





,19 


0, 


,046 


± 


0, 


,014 


0, 


,032 


± 


0, 


,014 


LP815-43 


2, 


,82 


± 


0, 


.25 


0, 


,036 


± 


0, 


,021 


0, 


,046 


± 


0, 


,022 
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Table 6. Predicted ^Li and ^Li depletion factors (in dex) for our 

PROGRAM stars INTERPOLATED FROM PREDICTIONS BY RICHARD ET AL. (2002, 

2005; AND 2005, private communication). The second and third columns 

GIVE THE EXPECTED DEPLETION DURING THE PRE-MAIN SEQUENCE (PMS) WHILE 
THE OTHERS ARE FOR THE EVOLUTION THEREAFTER TO WHICH THE PMS DE- 
PLETION FACTORS SHOULD BE ADDED. DIFFUSION DENOTES THE CASE OF PURE 
DIFFUSION WHILE T6.09 AND T6.25 ARE FOR TWO DIFFERENT MODELS OF TUR- 
BULENT DIFFUSION (see TEXT FOR DETAILS). 



Star 


PMS 


Diffusion 


T6.09 


T6.25 


m 


6Li 


^Li 


6Li 


^Li 


6Li 




6Li 


HD 3567 


0.01 


0.8 


0.3 


0.3 


0.2 


0.6 


0.5 


2.1 


HD 19445 


0.01 


0.9 


0.2 


0.2 


0.2 


0.6 


0.5 


2.2 


HD 59392 


0.02 


1.0 


0.2 


0.2 


0.2 


0.6 


0.5 


2.3 


HD 102200 


0.01 


0.8 


0.3 


0.3 


0.2 


0.6 


0.5 


2.0 


HD 106038 


0.02 


1.1 


0.2 


0.2 


0.2 


0.7 


0.6 


2.5 


HD 140283 


0.01 


0.4 


0.2 


0.3 


0.2 


0.4 


0.5 


1.9 


HD 160617 


0.01 


0.6 


0.4 


0.3 


0.2 


0.3 


0.5 


1.8 


HD 213657 


0.01 


0.5 


0.4 


0.4 


0.2 


0.4 


0.5 


1.8 


HD 298986 


0.01 


0.7 


0.3 


0.3 


0.2 


0.5 


0.5 


1.9 


HD 338529 


0.00 


0.3 


0.7 


0.7 


0.2 


0.4 


0.4 


1.7 


G 013-009 


0.00 


0.4 


0.6 


0.6 


0.2 


0.4 


0.4 


1.7 


G 020-024 


0.00 


0.4 


0.5 


0.5 


0.2 


0.4 


0.4 


1.8 


G 075-031 


0.01 


0.9 


0.3 


0.3 


0.2 


0.6 


0.5 


2.1 


G 126-062 


0.01 


0.5 


0.4 


0.4 


0.2 


0.4 


0.4 


1.8 


G 271-162 


0.00 


0.5 


0.5 


0.5 


0.2 


0.4 


0.4 


1.8 


BD+09°2190 


0.00 


0.3 


0.9 


0.8 


0.2 


0.3 


0.4 


1.6 


BD +03° 0740 


0.00 


0.4 


0.5 


0.5 


0.2 


0.4 


0.4 


1.7 


BD- 13° 3442 


0.00 


0.4 


0.6 


0.6 


0.2 


0.4 


0.4 


1.7 


CD -30° 18140 


0.00 


0.5 


0.5 


0.5 


0.2 


0.4 


0.4 


1.8 


CD-33°1173 


0.00 


0.3 


0.8 


0.9 


0.2 


0.3 


0.4 


1.6 


CD -33°3337 


0.02 


1.1 


0.2 


0.2 


0.2 


0.7 


0.6 


2.5 


CD -35° 14849 


0.00 


0.4 


0.5 


0.5 


0.2 


0.4 


0.4 


1.7 


CD -48° 2445 


0.00 


0.5 


0.5 


0.5 


0.2 


0.4 


0.4 


1.8 


LP 815-43 


0.00 


0.3 


0.9 


0.8 


0.2 


0.3 


0.4 


1.6 



